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ABSTRACT 

A  Computational  Study  of  the  Modification  of  Raindrop  Size 
Distributions  in  Subcloud  Downdrafts.  (August  1981) 

Robert  George  Borchers,  B.S.,  Texas  A&M  University; 

M.S.,  University  of  Michigan 

Co-Chairmen  of  Advisory  Committee:  Dr.  Phanindramohan  Das 

Dr.  James  R.  Scoggins 

A  computational  study  is  made  of  the  variation  of  steady-state 
raindrop-size  distributions  (0.004  cm  radius  <_  0.40  cm)  in  adiabatic 
subcloud  downdrafts  of  constant  magnitude.  The  cloud  base  drop-size 
distribution  is  assumed  to  be  Marshall-Palmer.  The  microphysical 
processes  progressively  introduced  are  evaporation,  collision- 
coalescence,  aerodynamic  and  collisional  breakup  of  drops.  Collision- 
coalescence  and  breakup  are  treated  through  a  stochastic  model. 
Thermodynamic  and  hydrometeoric  variables  of  temperature,  relative 
humidity,  total  liquid  water,  rainfall  rate  and  radar  reflectivity  are 
computed  as  functions  of  height  below  cloud  base  along  with  the  drop- 
size  distributions. 

It  is  found  that  evaporation  tends  to  deplete,  as  expected,  the 
smaller  members  of  the  droplet  population.  The  slope  of  the  Marshall- 
Palmer  distribution  as  determined  by  the  distribution  of  larger  size 
drops  remains  virtually  unaffected  by  evaporation.  When  collision- 
coalescence  is  included,  the  depletion  of  the  smaller  drops  is  enhanced 
by  the  larger  drops  sweeping  out  the  smaller  ones.  Aerodynamic  break¬ 
up  when  added  to  evaporation  and  collision-coalescence  has  little 
effect  on  the  resulting  drop-size  distribution  except  for  very  high 


initial  precipitation  rates.  When  collisional  breakup  is  added  to  the 
previously  mentioned  microphysical  processes,  the  effect  of  aerodynam¬ 
ic  breakup  is  virtually  eliminated  because  collisional  breakup  quickly 
depletes  the  concentration  of  drops  with  radii  greater  than  0.15  cm. 

Computed  raindrop-size  distributions  agree  quite  well  with 
measured  maritime  raindrop  spectra  reported  in  the  literature,  the 
agreement  being  better  for  the  higher  precipitation  rates  measured. 

Other  results  of  the  study,  as  anticipated,  are  (i)  that  the 
downdraft  remains  unsaturated  in  the  presence  of  precipitation,  the 
degree  of  subsaturation  increasing  with  the  strength  of  the  downdraft, 
and  (ii)  that  the  temperature  in  the  downdraft  lies  between  the  moist- 
adiabatic  and  the  dry-adiabatic  lapse  rates  of  temperature. 
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CHAPTER  I 


INTRODUCTION 

Raindrop  size  distribution  and  its  variation  in  the  vertical  is 
of  interest  because,  as  Srivastava  (1978)  pointed  out,  (i)  it  is 
important  in  understanding  processes  affecting  the  growth  and  evapora¬ 
tion  of  raindrops,  (ii)  a  knowledge  of  the  distribution  of  raindrop 
sizes  is  essential  for  remote  measurements  of  rainfall  by  radar,  and 
(iii)  the  size  distribution  of  raindrops  has  important  bearing  on  pre¬ 
cipitation  loading  and  its  effect  on  the  dynamics  of  clouds.  More 
specifically,  most  raindrop  measurements  are  made  at  the  ground  while 
the  precipitation  development  occurs  in  the  cloud.  Besides,  radar 
observations,  in  general,  relate  to  conditions  near  the  cloud  base.  It 
is  of  obvious  theoretical  and  practical  importance,  therefore,  to 
relate  the  drop-size  distribution  at  the  cloud  base  to  that  at  the 
ground . 

a.  Background 

The  modification  of  the  size  distribution  of  raindrops  during 
their  fall  has  been  investigated  by  several  researchers.  Early  in¬ 
vestigators  such  as  Mason  and  Ramanadham  (1954) ,  Rigby  et.  al.  (1954) , 
Sivaraman  and  Sivaramakrishnan  (1962) ,  and  Hardy  (1963)  calculated  the 
raindrop-size  distributions  resulting  from  mutual  coalescence, 
accretion  with  cloud  droplets,  and  evaporation  below  cloud  base. 

The  citations  on  this  and  the  following  pages  follow  the  style  of 
the  Journal  of  the  Atmospheric  Sciences. 
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They  assumed  an  exponential,  Marshall-Palmer  type  distribution  at 
the  melting  level,  no  downdraft,  and  a  constant  lapse  rate  of  tempera¬ 
ture  and  relative  humidity  below  cloud  base.  Their  results  indicated 
that  small  drops  were  depleted  by  each  modification  process  considered. 
The  concentration  of  large  drops,  on  the  other  hand,  increased  due  to 
collision-coalescence  and  accretion  but  decreased  by  evaporation. 

Hitschfeld  (1955)  and  Srivastava  (1967) ,  working  with  coalescence 
only,  concluded  that  an  initial  exponential  drop-size  distribution 
underwent  comparatively  little  change  between  cloud  base  and  ground. 

In  a  later  computation,  Srivastava  (1971)  added  aerodynamic  breakup 
of  water  drops  as  a  modification  process.  Using  the  data  of  Komabayasi 
et  al.  (1964)  on  the  probability  of  aerodynamic  breakup  and  the  size 
distribution  of  the  fragments  resulting  from  such  breakup,  he  found 
that  the  equilibrium  drop-size  distribution  was  similar  to,  but  flatter 
than,  observed  raindrop-size  distributions.  Srivastava  suggested 
that  the  effects  of  evaporation  and  collisional  breakup  may  help  pro¬ 
duce  better  agreement  between  theory  and  observation. 

Experiments  on  collisional  breakup  and  the  resulting  fragment- 
size  distributions  were  reported  by  Brazier-Smith  jet  al.  (1972) , 
McTaggart-Cowan  and  List  (1975)  and  Bradley  and  Stow  (1979) .  Brazier- 
Smith  et^  al.  (1973a,  b) ,  Young  (1975),  List  and  Gillespie  (1976), 
Srivastava  (1978)  and  Gillespie  and  List  (1978/79)  computed  the  evo¬ 
lution  of  drop-size  spectra  by  collisional  breakup  as  well  as  other 
processes.  Brazier-Smith  et^  al.  used  collision-coalescence  and 
collisional  breakup  to  compute  changes  in  drop-size  spectra.  They 
found  that  the  drop-size  distribution  after  several  minutes  was 
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insensitive  to  the  choice  of  an  initial  spectrum  and  became  markedly 
bimodal.  On  the  other  hand,  Gillespie  and  List  (1976)  considered  the 
same  two  modification  processes  and  found  that  the  initial  exponential 
drop-size  spectra  remained  essentially  exponential  as  the  drops  fell 
to  the  surface.  Gillespie  and  List  (1978/79)  extended  their  1976  work 
to  apply  to  a  vertical  rainshaft  in  a  steady  state  with  no  interac¬ 
tions  with  the  surroundings. 

Young  (1975)  and  more  recently  Srivastava  (1978)  considered  both 
colllsional  and  aerodynamic  breakup  together  with  collision-coalescence 
*nd  condensation  processes.  They  found  that  colllsional  breakup  dom¬ 
inated  over  aerodynamic  breakup.  The  resulting  equilibrium  raindrop 
distribution  was  in  good  agreement  with  the  empirical  Marshall-Palmer 
distribution.  Bradley  and  Stow  (1977)  developed  a  raindrop-size 
distribution  model  by  using  evaporation  and  drop  interaction  which 
incorporated  colllsional  breakup.  They  found  that  an  initial  exponen¬ 
tial  distribution  remained  essentially  exponential  between  cloud  base 
and  ground. 

Research  by  Kamburova  and  Ludlam  (1966)  showed  how  a  downdraft 
remained  unsaturated  in  spite  of  evaporation  of  the  rain  carried  with  it. 
Analogous  research  into  the  evaporation  of  raindrops  in  the  unsaturated 
downdraft  below  the  cloud  was  reported  by  Das  and  Subbarao  (1968,  1972), 
and  Caplan  (1969) .  They  assumed  the  precipitation  water  to  be  divided 
into  drops  of  uniform  size  falling  in  a  downdraft  of  constant  magnitude 
and  found  the  subsaturation  to  be  greater  with  smaller  water  content, 
larger  drop  sizes,  and  stronger  downdrafts.  The  temperature  in  the 
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downdraft  was  found  to  lie  between  the  values  obtained  by  either  a 
dry-adiabatic  or  moist-adiabatic  process.  Kintigh  and  Das  (1970) 
investigated  the  modifications  of  a  drop-size  distribution  as  it 
descended  from  cloud  base  to  ground  in  an  unsaturated  downdraft  by 
both  evaporation  and  collision-coalescence.  Comparing  their  evapora¬ 
tion  results  to  the  earlier  studies  of  monodisperse  precipitation, 
they  found  that  the  presence  of  smaller  drops  in  the  population  leads 
to  a  higher  relative  humidity  in  the  downdrafts. 


Obi ectives 


Except  for  the  incomplete  study  of  Kintigh  and  Das  (1970)  no 
effort  has  been  made  toward  incorporating  the  primary  physical  processes 
in  a  comprehensive  study  on  the  modification  of  drop-size  distribution 
in  a  downdraft.  The  objective  of  this  research  is  to  fill  this 
gap  by  constructing  a  raindrop-size  distribution  model  incorporating 
evaporation,  collision-coalescence,  aerodynamic  breakup,  and  collis- 
ional  breakup  of  drops  in  a  rainshaft  accompanying  a  steady-state, 
□onentraining  downdraft.  Calculations  will  be  made  for  five  dif¬ 
ferent  situations  in  order  to  isolate  the  modification  processes  which 
most  affect  the  raindrop-size  distribution.  These  situations  will  be 
characterized  by : 


a)  evaporation  only; 

b)  evaporation  and  collision-coalescence; 

c)  evaporation,  collision-coalescence,  and  aerodynamic  breakup; 

d)  evaporation,  collision-coalescence,  and  collisional  breakup; 


and 


e)  evaporation,  collision-coalescence,  collisional  breakup,  and 
aerodynamic  breakup. 

The  problem  as  posed  will  result  in  variable  lapse  rates  of  temperature 
and  relative  humidity  in  the  downdraft  based  on  the  amount  of  raindrop 
evaporation  in  the  subcloud  region.  Drop-size  distributions  for 
several  downdraft  velocities  and  precipitation  rates  will  be  calculated 
so  that  comparisons  with  observed  raindrop  spectra  can  be  made.  Ad¬ 
ditionally,  the  change  in  rainfall  rate  and  radar  reflectivity  will  be 
computed  in  the  subcloud  region.  An  effort  will  be  made  to  compare 
the  computational  results  of  this  research  to  measurements  of  raindrop- 
size  distributions. 


CHAPTER  II 


BASIC  EQUATIONS 

The  basic  equations  in  this  research  consist  of  a  thermodynamic 
equation,  a  moisture  equation  and  a  set  of  equations  for  drop  con¬ 
centration  including  the  effects  of  evaporation,  collision-coalescence 
and  aerodynamic  and  collisional  breakup  of  drops.  These  equations  are 
applied  to  the  subcloud  region  under  the  assumption  of  a  one¬ 
dimensional,  steady-state  downdraft. 


a.  Thermodynamic  Equation 

The  thermodynamic  equation  is  used  to  calculate  temperature  in 
the  downdraft  and  to  account  for  compressional  heating  as  well  as  the 
cooling  due  to  the  amount  of  liquid  water  evaporated  in  the  subcloud 
layer.  We  start  with  the  simplified  equation  derived  by  Das  (1969): 


d<fr 

dt 


-L  d£ 
_ v  _ 

T  dt 


* 


(1) 


Writing  d>  *  C^  In  9  +  constant,  we  have 

.  JEL  .  (2) 

dt  9  dt 

By  assuming  conditions  of  one-dimensionality  and  steady-state,  and 
employing  the  Poisson  relationship  for  potential  temperature,  the 
hydrostatic  equation,  and  the  equation  of  state,  we  can  write  the 
vertical  temperature  gradient  as 


*The  symbols  used  in  this  paper  are  listed  and  defined  in  Appendix  A. 


where  (=  g/C  )  is  the  dry-adiabatic  lapse  rate  of  temperature, 
b.  Equation  for  Water  Vapor 

The  second  term  on  the  right-hand  side  in  (3)  requires  the 
substantial  rate-of-change  of  water  vapor  mixing  ratio  with  time. 
The  equation  for  d£/dt  as  derived  by  Das  (1969)  is 


(4) 


where  N'  is  the  total  number  of  drops  per  unit  mass  of  air,  and  DM/Dt 
is  the  substantial  change  of  drop  mass  with  time  as  observed  by  an 
observer  moving  with  the  drop.  Note  that  (4)  assumes  a  monodisperse 
population.  This  substantial  change  is  explained  in  more  detail 
later.  For  a  polydisperse  population,  we  write  N^dr  for  the  number 
of  drops  per  unit  mass  of  air  lying  in  the  radius  interval  between 
radius  r  and  r  +  dr.  With  this  notation  Das  (1969)  has  generalized 
(4)  to 

dt  =  “  Ht  f  [Nrdr]  (5) 

Jo 


for  a  polydisperse  raindrop  population.  In  this  research  (5)  forms 
the  basis  for  computing  the  temporal  change  in  water  vapor  mixing  ratio 
needed  in  (3) . 
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c.  Equation  for  Drop  Concentration 

The  equation  for  drop  concentration  can  be  written  formally  as 


3n 


IT  +  V-  “r<V+V  ‘  <Vev  +  «r>CC  +  Wr>AB  +  <V 


CB 


(6) 


where  n^  is  the  number  of  drops  per  cubic  meter  per  radius  interval, 

V  is  the  velocity  of  the  gaseous  phase  (dry  air  plus  water  vapor) , 
is  the  velocity  of  the  liquid  phase  with  respect  to  the  gaseous  phase, 
and  Qr 's  represent  the  source  and  sink  effects  caused  by  the  respec¬ 
tive  microphysical  processes.  The  subscripts  EV,  CC,  AB,  and  CB 
respectively  represent  the  processes  of  evaporation,  collision- 
coalescence,  aerodynamic  breakup  and  collisional  breakup.  On  the 
left-hand  side  (LHS)  of  (6),  the  first  term  represents  the  local 
change  with  time  of  raindrop  concentration  and  the  second  term  is  the 
divergence  of  the  flux  of  raindrops. 

Using  the  continuity  equation,  dp /dt  +  pV-V  =  0,  and  noting  that 

N  =  n  /p ,  we  can  show  that 
r  r  ’ 

dN 

LHS  of  (6)  -  P~ +  V*(nrVr)  (7) 


and. 


v,(nrvr)  +  p  +  ^r^CC  +  ^r^AB  +  ^r^CB^  *  ^ 


Applying  the  assumption  of  one-dimensional,  steady-state  motion,  we 


arrive  at 
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CHAPTER  III 


MODIFYING  PROCESSES  FOR  RAINDROP-SIZE  DISTRIBUTION 

The  physical  processes  modifying  a  subcloud  drop-size  distribu¬ 
tion,  considered  in  this  research,  are  (a)  evaporation,  (b)  collision- 
coalescence,  (c)  aerodynamic  breakup,  and  (d)  collisional  breakup 
of  the  drops.  The  continuity  equation  for  raindrop  concentration  in 
(8)  symbolically  incorporates  these  processes.  An  equation  formu¬ 
lated  by  Das  (1969)  forms  the  basis  of  the  algorithm  used  in  this 
research  to  account  for  evaporation.  The  procedure  used  for  colli¬ 
sion-coalescence  is  similar  to  that  used  by  Ogura  and  Takahashi 
(1973).  The  treatment  of  aerodynamic  breakup  follows  the  develop¬ 
ment  of  Srivastava  (1971) ,  and  that  of  collisional  breakup  uses  the 
formulation  of  List  and  Gillespie  (1976).  The  details  of  these 
treatments  are  described  below. 

a.  Evaporation 


Considering  evaporation  to  be  the  only  modifying  process  in  a 
rainshaft  accompanying  a  steady,  one-dimensional  downdraft  of  magni¬ 
tude  wQ,  one  can  rewrite  (9)  as: 


dN 


dz 


E V 


-1 

wDP 


<s!,rY 


dz 


Wn  3r  V*  rdtj  • 


(12) 


In  this  expression,  the  first  term  on  the  right-hand  side  represents 
vertical  sorting  of  raindrops,  and  ~  is  the  growth  rate  of  the  drop 
due  to  evaporation  so  that  the  explicit  form  (3/3r)  (N^  dr/dt)  stands 


11 


for  the  symbolic  quantity  (Q_)-_,/p  in  (9).  In  other  words,  the  sec- 

r  ev 

ond  term  represents  the  decrease  in  drop  concentration  due  to  diffu¬ 
sion  of  vapor  away  from  the  drops  of  a  particular  radius  class. 

If  both  sides  of  (12)  are  multiplied  by  the  class  interval,  Ar, 
the  resulting  equation  can  be  expressed  as 


VT.r  dN  1  3_  (  dr\  Nr 
WD  dz  WD  3r  ^  dt/~  WD 


dr 

dt 


(13) 


where  N  *  N^Ar  and  p  and  V^,  are  assumed  to  be  constant.  Appendix 
B  provides  an  analysis  of  the  error  committed  in  assuming  constant 
p  and  VT-  The  error  incurred  does  not  significantly  affect  the 
results  since  we  look  at  the  relative  effects  of  each  microphysical 
process  on  the  drop-size  distribution.  Solving  for  the  change  of 
drop  concentration  with  height,  we  have 


(14) 


The  first  term  on  the  right  represents  the  shifting  of  the  drop 
categories  because  of  drop-size  changes  due  to  evaporation.  The 
second  term  represents  an  adjustment  due  to  change  in  the  radius 
intervals  for  the  categories. 

In  order  to  obtain  an  expression  for  the  rate  of  change  of  drop 
radius  given  in  (14),  a  spherical  raindrop  is  assumed.  The  mass  of 
a  given  drop  is  then 
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where  <5^  Is  the  density  of  liquid  water.  Taking  the  time  derivative 
of  both  sides  and  solving  for  the  rate-of-change  of  drop  radius  we 
have 


.  i  dM 

dr  _  _ r 

dt  "  4irri25  dt  *  (16) 

l 

dM 

As  pointed  out  by  Das  (1969),  — —  does  not  represent  all  the  mass 

dt 

change  of  a  drop  due  to  evaporation,  but  only  the  part  of  the  change 
moving  with  the  gaseous  parcel.  The  actual  change  of  mass  due  to 
evaporation  is  really  the  "substantial  change"  and  can  be  denoted  by 
the  "substantial  change"  and  can  be  denoted  by 


DM 
_ r 

Dt 


(17) 


Combining  (16)  and  (17),  we  have 

.  .  T  DM  dM  ~) 

H [dt  +  V  dr)  •  <18 

dM 

Now,  we  must  obtain  an  expression  for  so  that  (18)  can  be  used 
to  calculate  the  rate-of-change  of  drop  radius.  Imposing  the  condi¬ 
tion  of  a  one-dimensional,  steady-state  downdraft  on  (17)  one  can 
obtain 


dM 
_ r 

dz 


-1 

vv 


T.r 


DM 
_ r 

Dt 


(19) 


which,  when  combined  with  (18),  yields 
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dr 

dt 


(20) 


This  Is  Che  expression  used  to  compute  the  rate  of  change  of  drop 
radius  In  (14). 

DM 

The  mass  evaporation  rate,  — given  in  (20)  was  tabulated  by 
Kinzer  and  Gunn  (1951) .  Their  experimental  results  appear  in  the 
Smithsonian  Meteorological  Tables  (List,  1949)  for  raindrop  radii 
from  0.005  cm  to  0.21cm.  In  this  research,  the  raindrop  spectrum 
ranges  from  0.004  cm  to  0.4  cm,  which  requires  an  extension  of  the 
Kinzer  and  Gunn  results.  The  method  used  to  extend  the  Kinzer  and 
Gunn  raindrop  evaporation  data  is  shown  in  Appendix  C. 

According  to  Byers  (1965),  the  rate  of  change  of  mass  for  a 
spherical  drop  can  be  written  as 

HOl  =  ^r(S-_l) 

Dt  a '  +  b ' 


where 


R  T 
v 

e  D 
s 


and  S  is  the  saturation  ratio: 


S  -  e  /e 
a  s 


(22) 


Here  the  saturation  vapor  pressure,  eg,  is  computed  from  an  equation 


given  by  Murray  (1967): 


e  -  6.1078  exp  [17. 2694(T-273. 16) / (T-35 .86) ]  (23) 

s 

DM 

It  can  be  seen  from  (21)  that  is  negative  for  unsaturated  condi¬ 
tions  and  that  the  mass  evaporation  rate  increases  as  the  air  becomes 
more  unsaturated. 

b.  Collision-coalescence 


The  collision-coalescence  mechanism  is  commonly  named  the  collec' 
tlon  process  (Berry,  1967).  In  this  process  raindrops  grow  by  grav¬ 
itational  collection:  larger  drops  which  fall  faster  collide  and 
coalesce  with  smaller  drops  in  their  path.  The  importance  of  the 
collection  process  in  warm  rain  is  well  established.  Houghton  (1968) 
stated  that  even  under  conditions  in  which  the  ice  phase  is  respon¬ 
sible  for  the  initiation  of  precipitation,  collection  is  the  dominant 
growth  mechanism  in  the  nonfreezing  regions  of  the  cloud  and  is, 
therefore, important  in  shaping  the  drop-size  distribution. 

There  are  two  collection  models  used  in  precipitation  physics; 
they  are  commonly  referred  to  as  the  continuous  and  stochastic  models 
The  continuous  collection  model  is  a  straightforward  concept  and  can 
be  credited  to  such  early  workers  as  Schumann  (1939)  and  Findeisen 
(1939).  It  is  characterized  by  a  raindrop  distribution  consisting  of 
many  small  droplets  and  a  few  large  drops;  all  large  drops  grow  at  a 
continuous  rate  as  they  sweep  out  the  small  droplets.  The  basic  no¬ 
tion  underlying  the  stochastic  model,  on  the  other  hand,  is  that  of 
a  "fortunate  drop,"  and  was  ushered  in  by  Telford  (1955).  Berry 
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(1967)  described  the  essential  differences  between  the  continuous  and 
the  stochastic  models  and  his  interpretation  is  given  in  Fig.  1. 
According  to  the  stochastic  collection  model,  depending  on  an  assoc¬ 
iated  probability  of  collection,  a  fraction  of  the  large  drops  of  a 
given  class  will  collect  smaller  droplets  contained  in  the  swept-out 
volume  and  grow  into  a  larger  drop  class.  This  probabilistic  condi¬ 
tion  arises  from  the  random  locations  of  the  smaller  droplets  with 
respect  to  the  volume  swept  out  by  the  larger  drop.  During  any  parti¬ 
cular  time  interval,  only  a  fraction  of  the  drops  in  any  particular 
size  class  will  grow  into  larger  size  classes. 

The  stochastic  collection  model  is  the  one  employed  in  this  re¬ 
search.  Following  the  notation  of  Ogura  and  Takahashi  (1973),  the 
stochastic  collection  process  can  be  represented  by  the  symbolic 
integro-dlfferential  equation 

-  ijf  ¥w"''  1-t)dl  -  »i  £  » •  a.  t) «  «« 

where  i  and  !,  etc.,  refer  to  the  size  categories  of  the  droplets. 

The  first  term  on  the  RHS  of  (24)  represents  the  gain  in  due  to  drops 
of  radius  r^  coalescing  with  drops  of  radius  r^_^  to  form  drops  of 
radius  r^,  while  the  second  term  represents  the  loss  of  due  to  drops 
of  radius  r^  coalescing  with  all  other  drops. 

In  (24),  $(i,  Z)  is  the  collection  kernel  defined  in  terms  of  the 
3weep-out  rate  of  a  geometric  volume  relative  to  the  two  drops  of 
radius  r^  and  r^  and  the  collection  efficiency;  it  represents  the 
probability  that  two  drops  will  coalesce  per  unit  time.  For 


Fig.  1.  The  continuous  and  stochastic  models. 


A.  The  continuous  model  assumes  that  small  droplets  are 
swept  out  as  if  their  water  were  distributed  uniformly 
in  space,  and  that  100  large  drops  grow  at  the  same 
rate,  becoming  100  similar  larger  drops. 


B.  The  stochastic  model  assumes  that,  say  1  in  10  of  100 
large  droplets  will  collect  a  small  droplet  during  a 
given  time,  and  then  1  in  10  of  each  larger  size  will 
collect  a  smaller  droplet.  Large  droplets  then  grow 
at  different  rates  (after  Berry,  1967). 
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this  research,  the  collection  kernel  given  by  Srivastava  (1967)  is 
used ,  viz. , 

$(i,«.)  -  "(rj+r^)2  E(i,2.)|vi  -  vj  (25) 

where  r^  and  r^  are  the  radii  of  the  drops  with  mass  i  and  l,  respec¬ 
tively,  and  E(i,i)  is  the  collection  efficiency.  The  collection  effi¬ 
ciency,  E(i,£),  is  a  product  of  the  collision  efficiency  E^  (i,£)  and 
the  coalescence  efficiency  E2(i,i). 

In  this  study,  the  collision  efficiency  is  taken  as  equal  to  one 
representing  geometric  sweep-out.  In  the  experiments  revised  by 
List  and  Gillespie  (1976),  it  appears  that  for  raindrops  the  inertial 
effects  are  large  enough  to  prevent  the  incoming  small  drop  from  be¬ 
ing  deflected  noticeably  by  the  flow  around  the  larger  one.  Thus, 
raindrop  trajectories  are  essentially  straight  lines.  In  addition, 
in  the  portion  of  the  present  research  in  which  consideration  was 
given  to  collision-coalescence  in  conjunction  with  evaporation  and 
aerodynamic  breakup,  the  coalescence  efficiency  was  set  equal  to  one. 
This  was  done  so  results  could  be  compared  with  those  of  Srivastava 
(1971) .  When  collisional  breakup  was  included  with  collision- 
coalescence,  the  coalescence  efficiency  as  formulated  by  Whelpdale  and 
List  (1971)  was  used.  This  will  be  discussed  in  greater  detail 
later  in  this  section. 

An  important  constraint  on  the  stochastic  collection  process  is 
that  the  total  mass  under  distribution  should  remain  constant. 

This  serves  as  a  check  on  the  numerical  accuracy  of  the  finite  differ¬ 
ence  scheme.  This  mass  represents  the  liquid  water  content  of  the 
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rainshaft,  and  can  be  written  as 

=  f  nr  Mr  dr  (26) 

Jo 

c.  Aerodynamic  Breakup  of  Drops 

Langmuir  (1948)  proposed  a  chain  process  as  a  possible  mechanism 
for  producing  the  drop-size  distributions  observed  in  warm  clouds. 
According  to  Langmuir,  there  was  a  critical  drop  radius  beyond  which 
all  drops  would  Immediately  become  unstable  and  breakup.  He  adopted 
the  radius  of  0.3  cm  to  be  the  critical  value.  However,  the  experi¬ 
mental  work  of  Komabayasi  et  al.  (1964)  indicates  that  the  critical 
size  is  not  very  sharp;  instead,  a  gradual  growth  of  aerodynamic 
instability  appears  to  be  the  case.  Srivastava  (1971)  has  analyzed 
the  experimental  data  of  Komabayasi  et_  al.  (1964)  on  the  probability 
of  aerodynamic  breakup  of  water  drops  and  the  size  distribution  of 
the  fragments  resulting  from  the  breakup.  When  the  aerodynamic  break¬ 
up  is  the  only  process  modifying  the  drop-size  distribution,  the  rate 
of  change  of  with  time  can  be  represented  by 

(Si)AB-  (27> 

where  represents  the  probability  with  which  a  drop  of  radius  r^ 
breaks  up  per  unit  time,  and  Q(i,Jl)  represents  the  number  of  frag¬ 
ments  of  radius  r^  given  by  the  breakup  of  one  drop  of  radius  r^. 
Expressions  for  and  Q(i,£)  are 


19 


(28) 


(29) 


where  is  expressed  in  centimeters.  The  expressions  for  and 
Q  (i,  £)  as  given  above  insure  the  requirement  for  total  mass  conser¬ 
vation. 


Collisional  Breakup  of  Drops 


On  reviewing  the  experimental  studies  on  aerodynamic  breakup 
including  those  of  Komabayasi  et.  al.  (1964) ,  List  et.  al.  (1970)  im¬ 
plied  that  aerodynamic  breakup  did  not  appear  to  explain  the  ob¬ 
served  size  distributions  of  raindrops.  Computations  of  Srivastava 
(1971),  who  regarded  aerodynamic  breakup  to  be  the  only  drop- 
disintegration  process,  produced  a  raindrop  size  distribution 
having  relatively  more  large  drops  than  observed  in  rain  formed  in 
all-liquid  processes.  Consequently,  in  the  Langmuir  chain  model 
of  warm  rain  development  a  more  effective  drop-disruption 
process  appears  to  be  needed  and,  according  to  List  et  al.  (1970), 
the  process  of  temporary  collision,  between  drops  of  sizes 
too  small  to  be  aerodynamically  unstable,  appears  to  be  the  answer. 
This  process,  commonly  referred  to  as  collisional  breakup,  has  been 
studied  experimentally  by  List  et_  al.  (1970),  Brazier-Smith  et  al. 
(1972),  McTaggart-Cowan  and  List  (1975),  and  Bradley  and  Stow  (1979). 
List  et  al.  (1970)  investigated  the  fragment  size  spectra  developing 


after  collisions  of  drops  with  radii  0.10  to  0.23  cm,  but  did  not 
consider  drops  falling  at  terminal  velocity.  They  concluded  that,  in 
order  to  study  the  growth  of  warm  rain  by  chain  process  and  to  make 
definite  statements  about  drop-size  limitations,  more  collision  and 
breakup  experiments  would  be  needed  in  a  larger  range  of  smaller  and 
larger  drop  combinations,  preferably  with  drops  falling  at  terminal 
velocity.  Experimental  studies  by  Brazier-Smith  et_  al.  (1972)  in¬ 
volving  drops  of  radii  0.015  to  0.075  em  and  later  studies  by  Brad¬ 
ley  and  Stow  (1979)  involving  drops  in  the  radius  range  0.06-0.17 
cm  failed  to  use  drops  falling  at  terminal  velocity.  McTaggart- 
Cowan  and  List  (1975)  updated  the  studies  of  List  et  al.  (1970)  by 
conducting  experiments  with  drops  of  radii  0.05  to  0.23  cm,  all 
falling  at  terminal  velocity.  They  used  a  high-speed  photographic 
method  to  observe  several  hundred  collisions  of  five  different  drop 
pairs.  The  radii  (rg)  of  the  smaller  drops  were  0.05  or  0.09  cm 
while  those  of  the  large  ones  (r^)  were  0.15,  0.18,  or  0.23  cm. 

Their  work  showed  that  the  coalescence  of  large  drops  with  smaller 
ones  was  a  rare  event,  the  most  likely  outcome  of  the  collisions  be¬ 
ing  the  production  of  numerous  fragments  with  loss  of  mass  from  the 
large  participating  drop.  Thus  the  coalescence  efficiency  was  much 
less  than  unity. 

Since  McTaggart-Cowan  and  List  (1975)  experimented  with  drops 
at  their  terminal  velocities,  it  is  their  fragment  size  distribu¬ 
tions  that  have  been  chosen  for  this  research.  Using  their  experi¬ 
mental  observations,  List  and  Gillespie  (1976)  developed  formal 
expressions  for  the  statistical  rates  of  production  of  small  and 


21 


large  fragments  resulting  from  collisional  breakup  of  a  pair  of 
drops.  The  contribution  of  these  fragments  to  a  drop-size  category 
is  expressed  through  a  function  W^(L,s)  which  is  the  rate  of  formation 
of  drops  of  radius  r^  by  collision-induced  breakup  of  a  larger  drop 
of  radius  r  interacting  with  smaller  drops  of  radius  rs: 

L 

/*k  2 

W  (L,s)-pirl  Ns(rL+rs)  E1(L,s)(VL-Vg)Pi(L,s)ds  .  (30) 

Jq 

Here  P.(L,s)  denotes  the  probability  of  obtaining  a  drop  of  radius 
r^  by  collision  between  and  subsequent  breakup  of  a  larger  drop  of 
radius  rT  and  a  smaller  drop  of  radius  r  .  Thus  the  rate  of  increase 
of  drop  concentration  in  i-th  category  can  be  written 


(w  +V 
^  D  T 


N.W.dr 
L  i  L 


“idrL 


(31) 


where  is  the  drop  concentration  per  unit  mass  of  air.  The  first 
term  on  the  R.H.S.  of  (31)  represents  a  gain  in  the  i-th  cate¬ 
gory  due  to  collision  between  drops  of  radii  larger  than  r^,  whereas 
the  second  team  gives  the  loss  in  the  category  due  to  their  collison 
with  smaller  drops. 

In  the  experiments  of  McTaggart-Cowan  and  List*  a  typical  over¬ 
all  distribution  of  fragments  shows  two  peaks,  one  toward  the  small 
and  the  other  toward  the  large  drop  end.  List  and  Gillespie  conclude 
that  the  observed  fragment  distributions  are  fairly  well  represented 
by  the  following  expressions: 
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f  =  1.44  x  103  (r3r  )1/2  (0.41  -  0.30  —  )  ,  (32) 

L  s  rL 


D  -  -1.0  - 


0.0196 


(33) 


with 


P(s)  =  (0.3^)  (1.529  x  105M1)° 


.  if  M  <_  1  x  10  6Kg,  (34) 


or 


7.  (M  x  106) 

P(s)  =  0.3  f )  — - r- 

rs  (0.0654)® 


,  if  M  >  1  x  10"6  Kg,  (35) 


where  f  is  the  total  number  of  small  fragments,  and  D  is  the  power 
law  constant  used  to  fit  the  small-fragment  distribution.  In  these 
expressions  the  radius,  r,  and  mass,  M,  are  expressed  in  units  of  centi¬ 
meters  and  kilograms,  respectively. 

A  Gaussian  function  fitted  by  List  and  Gillespie  to  the  large 
fragment  distribution  is  written* 


fH  exp  C-(Mi-Ml)2H2/2] 


P(L)= 


for  <_  (^  +  Mg) 
for  M±  >  (^  +  Ms) 


(36) 


with 


*There  is  a  sign  error  in  the  expression  as  printed  and  it  is 
possible  that  the  computation  of  List  and  Gillespie  (1976) 
failed  to  correct  the  expression;  further  comments  on  this  possi¬ 
bility  are  made  later  in  this  dissertation. 
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0.37 _ _ 

r  4r  (0.41-0.3r  /r  ) 

S  L  S  L 


(37) 


The  distribution  is  cut  off  for  +  Mg  because  no  fragment  lar¬ 
ger  than  the  sum  of  the  incoming  masses  can  be  produced.  Finally,  the 
overall  fragment  probability,  P  (L,s)  is  written 

Pi(L,s)  *  P (L)  +  P(s)  •  (38) 


As  mentioned  earlier,  the  experiments  of  McTaggart-Cowan  and 
List  (1975)  indicated  the  coalescence  efficiency  of  interacting  drops 
to  be  much  less  than  unity.  For  the  purpose  of  this  research,  the  coal¬ 
escence  efficiency  according  to  Whelpdale  and  List  (1971)  is  used, 
modified  to  agree  with  the  absence  of  coalescence  once  the  small 
drop  exceeds  0.05  cm  radius.  The  resulting  form  of  the  coalescence 
efficiency,  E2>  is 


E2  (L,s) 


for  rg  <_  0.05  cm 

for  r  >  0.05  cm 
s 


(39) 


It  is  understood  that  such  clearly  defined  limits  as  those  in 
(39)  probably  do  not  occur  in  nature.  These  coalescence  efficiencies 
are  used  so  our  results  can  be  compared  with  those  of  List  and  Gilles¬ 
pie  (1976)  and  Gillespie  and  List  (1978/1979). 


CHAPTER  IV 


PROCEDURE 

The  change  in  raindrop  concentration  in  a  subcloud  downdraft 
results  from  the  microphysical  processes  incorporated  in  (14),  (24), 
(27),  and  (31).  This  study  consists  of  a  numerical  integration  of 
these  four  equations  as  an  initial  value  problem.  Conditions  at  the 
cloud  base  consist  of  saturated  air  at  an  assumed  temperature,  a 
Marshall-Palmer  drop-size  distribution  corresponding  to  an  assigned 
rainfall  rate,  and  a  constant  subcloud  downdraft.  As  for  the 
choice  of  the  Marshall-Palmer  distribution  at  the  cloud  base, we  may 
recall  that  this  expression  provided  a  good  average  fit  to  the  exper¬ 
imental  data  of  Laws  and  Parsons  (1943)  except  for  a  tendency  to 
overestimate  the  numbers  of  small  droplets  (r  <  0.001  cm).  Although 
the  measurements  of  Laws  and  Parsons  were  made  at  the  surface,  the  work 
of  several  researchers  suggests  that  an  exponential  drop-size  dis¬ 
tribution  may  be  valid  aloft.  Srivastava  (1967),  who  summarized 
previous  drop-size  distribution  studies,  suggests  that  the  exponential 
distribution  may  be  quasi-stable  with  respect  to  modification  during 
fall,  and  while  the  exponential  distribution  itself  changes  slowly, 
other  distributions  might  tend  toward  an  exponential  shape.  In 
addition,  we  want  to  compare  our  computed  drop-size  distributions  with 
those  obtained  in  other  studies  that  used  the  Marshall-Palmer 
distribution  at  cloud  base. 

In  the  computational  setup,  the  subcloud  region  is  divided  into  a 
finite  number  of  layers  and  the  thermodynamic  and  hydrometeoric 


quantities  are  calculated  at  the  base  of  eacy  layer  for  a  given  subcloud 
downdraft  and  a  set  of  subcloud  conditions.  The  computation  is  con¬ 
tinued  down  to  1500  m  below  the  cloud,  which  is  referred  to  as  the 
"surface."  Once  the  temperature,  relative  humidity,  total  liquid 
water,  rainfall  rate,  radar  reflectivity,  and  drop-size  distribution 
at  the  surface  are  calculated  for  a  given  set  of  conditions,  the  pro¬ 
cedure  is  repeated  for  other  constant  downdraft  velocities  and  cloud 
base  precipitation  rates  in  order  to  gauge  the  sensitivity  of  the 
resulting  steady-state  distributions  to  these  parameters.  Our 
study  uses  constant  downdraft  velocities  of  5,  10,  and  15  m/sec 
with  cloud-base  precipitation  rates  of  25,  50,  75,  and  100  mm/hr. 

a.  Cloud-Base  Conditions 

The  Marshall-Palmer  distribution  of  raindrops  at  the  cloud  base 
is  written 

n^  =  0.16  exp  (-Ar),  A  =  82  R  0.21  (40) 

-4 

where  nf  has  units  of  cm  ,  r  is  in  cm,  and  R  is  in  mm/hr.  Fig. 

2  shows  the  Marshall-Palmer  drop-size  distributions  for  several 
rainfall  rates. 

For  numerical  purposes  the  drop-size  distributions  are  discre¬ 
tized  into  40  categories  with  mean  radii  (r^  i  *  1,40)  ranging  from 
0.004  to  0.4  cm  on  a  logarithmic  scale.  Table  1  gives  n^  according  to 
(40)  for  each  radius  category  along  with  the  mass  of  each  drop  in  that 
category  as  well  as  the  radius  interval  (Ar^)  for  each  category.  The 
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Table  1.  Marshall-Palmer  drop-size  distribution  for  a  precipi¬ 
tation  rate  of  50  nun/hr. 


r. 

l 

Vi 

N. 

i 

M. 

l 

(cm) 

(cm) 

(m/sec) 

(/// m^cm) 

(kg) 

.0042 

.0005 

0.24 

1.373x10? 

3.187x10“}° 

4.507x10“}!} 

6.374x10“}° 

.0048 

.0006 

0.26 

1. 348x10? 

.0053 

.0006 

0.29 

1.320x101 

.0060 

.0007 

0.33 

1. 289x10c 

9.014x10“?° 

.0067 

.0008 

0.38 

1. 255x10c 

1.275x10“? 

.0076 

.0009 

0.44 

1.218x10? 

1.803x10“: 

.0085 

.0010 

0.52 

1.178x10? 

2.549x10“: 

.0095 

.0011 

0.63 

1.135x10* 

3.605x10? 

.0109 

.0012 

0.73 

1.088x10* 

5.099x10“? 

.0120 

.0014 

0.84 

1.038x10* 

7.211x10* 
1.020x10“: 
1.422x10“: 
2.039x10“® 
2.884x10“: 
4.079x10“: 
5.768x10 “J 
8.157x10, 

.0135 

.0016 

0.97 

9.847x10* 

.0151 

.0018 

1.10 

9.279x10* 

.0170 

.0020 

1.21 

8.680x10* 

.0190 

.0022 

1.43 

8.054x10* 

.0214 

.0025 

1.62 

7.404x10* 

.0240 

.0028 

1.84 

6.738x10* 

.0269 

.0031 

2.08 

6.060x10. 

5.381x10* 

4.709x10* 

.0302 

.0035 

2.33 

1.154x10", 

,0339 

.0039 

2.63 

1.631x10, 

.0381 

.0044 

2.93 

4.054x10* 

2.307x10", 

.0427 

.0049 

3.29 

3.427x10* 

3.263x10, 

,0479 

.0055 

3.66 

2.837x10* 

4.614x10", 

.0538 

.0062 

4.04 

2.296x10* 

6.525x10, 

,0604 

.0070 

4.43 

1.810x10* 

9.228x10"? 

,0678 

.0078 

4.82 

1.386x10* 

1.305x10"? 

,0761 

.0088 

5.22 

1.027x101 

1.846x10"? 

0854 

.0099 

5.66 

7.339x10* 

2.610x10"? 

,0959 

.0111 

6.07 

5.032x10* 

3.691x10"? 

,1076 

.0124 

6.53 

3.294x10* 

5.220x10"? 

,1208 

.0140 

7.00 

2.048x10* 

7.382x10"? 

,1356 

.0157 

7.43 

1.201x10* 

1.044x10"* 

1522 

.0176 

7.87 

6.596x10, 

1.476x10"? 

1708 

.0198 

8.23 

3.367xl0£ 

2.088x10"* 

1918 

.0222 

8.58 

1.583xl02 
6.784x1  o} 

2.953x10'® 

,2152 

.0249 

8.79 

4.176x10"? 

2416 

.0279 

9.00 

2.621x10? 

5.905x10"* 

2712 

.0314 

9.08 

9.013x10° 

8.351x10": 

3044 

.0352 

9.16 

2.720x10° 

1.181x10"* 

3417 

.0395 

9.20 

7.087x10“ 1 

1.670x10'* 

3835 

.0443 

9.25 

1.566x10"* 

2.362x10'* 

4305 

.0498 

9.27 

2.877x10“^ 

3.340x10"* 

method  used  to  determine  the  radius  intervals  of  the  categories  is  given 
later.  Fig.  3  shows  the  percent  of  total  liquid  water  contained  in  each 
drop  radius  category  for  the  Marshall-Palmer  distribution  at  a  rainfall 
rate  of  50  mm/hr.  It  is  easily  seen  that  the  bulk  of  liquid  water  is 
contained  in  the  drop  radii  ranging  from  0.004  to  0.4  cm. 

A  temperature  of  5°C  is  assumed  at  cloud  base.  As  already 
mentioned,  constant  downdrafts  of  5,  10,  and  15  m/sec  are  specified 
at  the  cloud  base  and  kept  constant  in  the  subcloud  layer. 

b.  Grid  System  and  Numerical  Procedures 

The  subcloud  region  is  assumed  to  be  1500  m  deep  and  is  divided 
into  layers  of  25-m  thickness.  Consequently 

z  -  JAz,  (41) 

where  j  *0,  1,  2,  . ..,  61  and  Az  ■  25  m,  with  j  *  0  at  the  surface 
and  j  ■  61  at  the  cloud  base. 


c.  Numerical  Formulation  for  Raindrop  Size 

The  size  spectrum  of  raindrops  is  divided  into  40  logarithmic 
categories  covering  the  radius  range  0.004-0.4  cm  in  a  manner  similar 
to  Jorgensen  (1974)  with  droplets  of  r  <  0.004  cm  lumped  into  one 
group.  The  conversion  to  the  logarithmic  scale  is  made  so  that  the 
number  of  data  points  would  be  computationally  manageable  and,  as 
will  be  seen  later,  the  calculation  of  raindrop  mass  would  be  simpli¬ 
fied.  In  order  to  facilitate  the  use  of  the  forward  difference  scheme 
for  the  40th  drop  category,  one  extra  category  is  added  at  the  larger 


raindrop  distribution  for  a  precipitation  rate  of  50  mm/hr. 


f 
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drop-size  end  of  the  distribution.  Starting  with  the  smallest  value 
of  0.004  cm,  the  drops  fall  into  size  classes  that  are  obtained  from 
the  relation 


k+1 


10(lo810rk 


+  DLOGR) 


1.  2, 


40 


(42) 


where  DLOGR  is  defined  as  the  difference  between  the  logarithms  of 
radii  of  any  two  adjacent  size  classes  and  is  written 

DLOGR  =  (log102)/6  . 

Berry  (1967)  introduced  the  logarithmic  radius  intervals  as  used 
in  this  study  in  order  to  emphasize  the  resolution  in  the  smaller 
size  classes  which  can  be  shown  as  follows.  The  mass  of  a  drop  in 
the  k-th  class  is 


Mk  ■  T  ■  (43> 

The  mass  of  the  (k+l)-th  class  can  be  shown  to  be 

"kH  "  Mk  103  DL0GR  -  “fc-l  ’ 

which  means  that  the  mass  of  the  raindrops  doubles  every  two  radius 
classes.  The  high  resolution  at  the  smaller  size  end  of  the  spectrum 
is  especially  suitable  for  this  research  since  the  small  particles 
change  size  rapidly  under  the  effect  of  evaporation. 

The  drop  radius  calculated  from  (42)  represents  the  midpoint 
radius  of  a  particular  drop  category.  Table  1  gives  the  midpoint  ra¬ 


dii  r^,  and  the  radius  intervals,  Ar^  (in  cmX  over  which  these  midpoint 


31 


32 


C  1+1  P1+l 
ea,j+l  “  0.622  +  £ 


(48) 


where  P  is  the  pressure  at  the  new  level.  The  pressure  is  computed 


from  the  hydrostatic  relationship 


P  -  P 
j+1 


fgAz 

i  ei  v 

— 


(49) 


where  T  is  the  mean  temperature  in  the  layer. 


e.  Numerical  Formulation  of  the  Evaporation  Effects 


In  formulating  the  finite-difference  approximation  to  (14) 

we  proceed  as  follows.  Writing  N  “  (N  Ar).  where  N,  represents  the 

k.  r  k  k 

concentration  of  drops  in  the  category  k,  we  can  rewrite  (14)  as 


_3_ 

3r 


(Nr3T4')  '  »r  jt  h  (4r>] 


(50) 


in  which  the  last  term  arises  from  the  fact  that  the  categories 
have  variable  width  Ar^.  The  finite-difference  form  of  (50)  with 
respect  to  r  can  be  written 


^k  (l _ \ 

dz  EV  \wD+VTj 


(ArNr,dr/dt)k+1  -  (ArNrdr/dt)k 
rk+l  "  rk 


k+1  ~  Ark  1 

k+1  '  rk  _J 


We  observe  the  following: 
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Ark  "  rk+l  -  rk  =  rk(Y  “l)* 


Ark+1  “  rk+2  '  rk+l  ”  rk+l 


where 


rk+2/rk+l  ”  rk+l/'rk 


which,  because  of  constant  logarithmic  intervals  in  our  categoriza¬ 
tion,  is  a  constant.  Using  these  results  we  can  reduce  (51)  to 

®l  *L  ■  (■■  »), J  - 

Rewriting  (52)  in  terms  of  finite-differences  in  z,  we  have 

W  [k  - 


(\,j+i)ev  ■  "k,J  *  <4V 


EV 


Eqs.  (53)  and  (54)  are  used  to  compute  the  drop-size  distribution 
which  reflect  the  change  due  to  evaporation  only.  Once  the  number 
concentration  in  each  radius  category  is  computed  ^t  the  new 
level  j+1,  the  liquid-water  content  at  that  level  may  be  computed 


\}+l  "kJ1  MkNk,j+lp 


where  is  mass  of  one  spherical  drop  in  the  radius  category  k. 
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f.  Numerical  Formulation  for  Collision-Coalescence 


A  formulation  similar  to  that  given  by  Ogura  and  Takahashi 
(1973),  is  used  to  change  the  collision-coalescence  equation  given  in 
(24)  to  a  summation  suitable  for  programming: 


wVcc 


r  3 

Ip  Az  \  kh  /r,  \ 

Ivv  k,E=1  v Nkc  (^r) 


kmax 

Z 

k'  -1 


4>(kc,k?)  -  N,  Z  N  *(k,k’  ) 


or 


(56) 


5  <4VcC-fl  -  «\>CC-L 


Here  the  subscripts  CC-G  and  CC-L  stand  for  gain  and  loss  respective¬ 
ly,  in  the  concentration  of  drops  in  the  k  category  due  to  collision- 
coalescence.  In  addition, 

k 

kh  =  k - —  In  2  =  k  -2  k  DLOGR  In  10  ,  (57) 

3  o 

where  kQ  is  an  adjustable  scale  factor  that  fixes  the  desired  data 
points  at  integer  values  of  k.  Since  kQ  =  1/ (DLOGR  In  10),  kh  ■  k-2, 
which  is  the  class  index  of  drops  whose  mass  is  half  that  of  class 
k.  Further 

kc  "  k  +  T  ln  jj-  "  exp  3(k'-k)/kJ  .  (58) 

To  calculate  N^,  the  following  four-point  Lagrangian  interpolation 


scheme  is  used: 
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Nkc  "  “  i  (l-a)(a)(l+a)Nk_2  +  j  (2-a)  (a)  (1+cON^ 
+  j  (2-a) (1-a) (l+a)Nk  -  \  (2-o) (1-a) (a)Nfc+1 


(59) 


where 


a 


k-k 

c 


L  3k'  -k  1 

L  p  ~  J 


(60) 


Once  (ANk)cc  is  calculated,  the  number  concentration  in  the  k- 
category  is  obtained  by  combining  the  effect  of  evaporation  and 
collision-coalescence,  that  is. 


^k.j+l^EV+CC 


Nk,j  +  (ANk}EV  +  (ANk>CC 


(61) 


The  collision-coalescence  process  operates  under  the  constraint 
of  conservation  of  total  liquid  water.  No  drops  are  allowed  to 
enter  the  drop-size  distribution  at  the  small  radius  end  (<  0.004  cm) 
due  to  collision-coalescence,  but  drops  are  allowed  to  move  out  the 
larger  radius  end  (>  0.4  cm)  of  the  spectrum.  The  first  term  in 
brackets  in  (56)  is  allowed  to  operate  for  the  41st  radius 
category.  This  allows  drops  to  collect  in  th^-s  category  due  to  col¬ 
lection  in  the  smaller  radius  categories,  but  does  not  allow  any 
drop  to  leave  due  to  collection  between  category  41  drops  and  all 
other  drops  in  the  distribution.  Although  this  does  not  conserve 
total  liquid  mass  between  drops  of  radii  0.004  cm  to  0.40  cm,  it  does 
conserve  mass  when  the  aerodynamic  breakup  process  is  considered. 


g.  Numerical  Formulation  for  Aerodynamic  Breakup 

The  finite  difference  formulation  of  Jorgensen  (1974)  is  used 
to  evaluate  the  aerodynamic  breakup  term  given  in  (27),  i.e.. 


<4Vab  -  (  » 


Az 


r 


D+VT,k< 
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-NkPk+  E 


k  -k+1 


N  Q(k,k’)P 
k  k' 


(62) 


where  the  functions  P  and  Q  are  given  by  (28)  and  (29),  respectively. 
On  combining  the  effects  of  evaporation,  collision-coalescence,  and 
aerodynamic  breakup,  we  have 


‘Vj+lW+AB  -  \,j  +  (SNk>EV  +  (4VcC  +  <4VaB  •  (63) 


h.  Numerical  Formulation  for  Collisional  Breakup 

The  numerical  procedure  of  List  and  Gillespie  (1976)  and  Taka- 
hashi  (1977)  was  reviewed  in  formulating  the  collisional  breakup 
procedure.  The  collisional  breakup  computation  is  divided  into  two 
parts;  one  for  the  gain  term  (subscripted  CB-G)  and  the  other  for  the 
loss  term  (subscripted  CB-L).  The  gain  term  is  computed  from 


^ANk^CB-G 


PaAz 


41 

£ 


L-l 

l 


(wD+VT,k)  1=32  s=23 


W(W  VL’S>  * 


(64) 


Eq.  (64)  represents  the  gain  in  the  concentration  of  drops  in  the 
radius  category  k  due  to  breakup  following  collision  between  drops 
in  the  radius  category  r  (large)  and  r  (small).  As  mentioned 

Li  S 

previously,  the  smallest  r^  used  by  McTaggart-Cowan  and  List  (1975) 
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was  0.15  cm  while  the  smallest  r  was  0.05  cm.  The  minimum  rT  and 

s  L 

rs  correspond  to  approximately  categories  32  and  23,  respectively. 
Although  their  experimental  results  did  not  cover  the  entire  drop 
size  distribution  of  our  research,  the  P(L)  and  P(s)  functions  were 
considered  continuous  since  they  approached  zero  at  r  *  r  y  and 

r  =  r  .  .  The  one  exception  for  P(s)  was  when  the  mass,  M,  ,  was  less 
turn  k. 

than  1  x  10  ^  kg.  Rather  than  let  P(s)  increase  exponentially  at 

the  small  radius  end  of  the  spectrum,  it  was  considered  constant  at 

r^  =»  0.025  cm  (radius  category  16).  This  was  done  because  McTaggart- 

Cowan  and  List  (1975)  could  not  resolve  drop  fragment  information 

below  this  radius.  Their  results,  however,  appear  to  be  reaching 

a  fragment  maximum  near  this  radius  value. 

The  loss  of  drops  in  the  radius  category  k  is  considered  in  three 

different  ways  depending  on  whether  (i)  k  <  23,  (ii)  23  <_  k  <  32,  or 

(iii)  k  ^  32.  No  collisional  breakup  is  allowed  when  k  <  23 

(r^  <  0.05  cm).  This  is  consistent  with  the  experimental  data  of 

McTaggart-Cowan  and  List.  When  23  <_  k  <  32  (0.05  cm  <  rfc  <  0.15  cm), 

the  loss  term  for  the  drops  in  category  k  is  expressed  as 

<iNk,CB-L'( 

(65) 

where  fr  refers  to  a  fragment  drop  resulting  from  the  collisional 
breakup  of  drops  of  radius  r^  and  r^.  In  order  to  determine 
the  number  of  drops  in  category  k  lost  by  collision  with  category 
L  drops,  the  total  mass  of  fragments  produced  is  multiplied 


P  A z  .  41  N^  40 

:r+v — )  1  m~st  s  $(rT.r1)pfra,k) 

WD+VT,k'  L-32  fr-1  fr  k  L  k  fr 
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The  numerical  factors  in  the  above  expressions  appear  because  of 
conversion  to  the  conventional  units  from  the  units  used  in  the  com¬ 
putation,  namely,  VT  ^  in  m/sec,  and  r^  in  cm.  Also  note  that  the 
conventional  Z  is  computed  with  diameters  rather  than  radii  of  drops. 
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CHAPTER  V 

RESULTS 


In  presenting  the  results  we  recall  our  objective  which  is  to 
isolate  the  effects  of  the  microphysical  modification  processes  of 
evaporation,  collision-coalescence,  aerodynamic  breakup,  and  collis- 
ional  breakup  on  the  drop-size  distribution  in  a  rainshaft  accompany¬ 
ing  a  steady-state,  non-entraining  downdraft.  To  that  end  computations 
were  made  for  five  cases  as  follows: 

Case  A:  Evaporation; 

Case  B:  Evaporation  and  collision-coalescence; 

Case  C:.  Evaporation,  collision-coalescence,  and  aerodynamic 
breakup ; 

Case  D:  Evaporation,  collisional-coalescence,  and  collisional 
breakup ;  and 

Case  E:  Combination  of  all  the  processes  mentioned  in  Case  A 
through  Case  D . 

The  variations  of  drop-size  distribution  resulting  in  each  case  are 
presented  along  with  variations  in  temperature,  relative  humidity, 
rainfall  rate,  and  radar  reflectivity  in  the  subcloud  region.  Compu¬ 
tations  were  made  by  using  steady  state  downdrafts  of  5,  10,  and  15  m/ 
3ec,  and  for  still  air  precipitation  rates  of  25,  50,  75  and  100  mm/hr. 


The  vertical  distributions  of  selected  thermodynamic  and  moisture 
variables  resulting  from  evaporation  are  listed  in  Tables  3,  4,  and  5 
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Table  3.  Vertical  distributions  of  thermodynamic  and  hydrometeoric 
quantities  for  Case  A  in  a  constant  subcloud  downdraft  of 
5  m/sec,  where  z  ■  height  above  the  surface  (1500  m  ■  cloud 
base),  TD  «  dry-adiabatic  temperature,  T  »  rainshaft  temp¬ 
erature,  S  ■  relative  humidity,  M  ■  liquid  water  content, 

R  ■  rainfall  rate,  and  Z  *  radar  reflectivity. 


z 

(m) 

td 

(K) 

T 

00 

S 

(7.) 

MCxlO-5) 

(kg/m3) 

R  Z(xl0‘5) 

(mo/hr)  (n*n6/m3) 

R  =* 

25  mm/hr 

1500 

278.0 

278.0 

100.0 

1.24 

25.0 

0.33 

1000 

282.9 

282.8 

77.4 

1.22 

24.7 

0.33 

500 

287.8 

287.3 

62.4 

1.17 

24.0 

0.32 

0 

292.7 

291.6 

51.8 

1.10 

22.8 

0.31 

R  =» 

50  mm/ hr 

1500 

278.0 

278.0 

100.0 

2.24 

50.0 

0.92 

1000 

282.9 

282.7 

78.1 

2.21 

49.7 

0.91 

500 

287.8 

287.0 

64.3 

2.15 

48.6 

0.90 

0 

292.7 

291.1 

54.9 

2.06 

46.9 

0.89 

R  = 

75  mm/hr 

1500 

278.0 

278.0 

100.0 

3.28 

75.0 

1.77 

1000 

282.9 

282.6 

78.6 

3.25 

74.5 

1.76 

500 

287.8 

290.7 

65.9 

3.17 

73.1 

1.74 

0 

292.7 

290.7 

57.6 

3.06 

71.1 

1.71 

R  = 

=  100  mm/hr 

1500 

278.0 

278.0 

100.0 

4.03 

100.0 

2.51 

1000 

282.9 

282.6 

79.0 

3.99 

99.4 

2.50 

500 

287.8 

286.7 

66.9 

3.90 

97.9 

2.47 

0 

292.7 

290.5 

59.3 

3.79 

95.6 

2.44 
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Table  4.  Vertical  distributions  of  thermodynamic  and  hydro- 
meteoric  quantities  for  Case  A  in  a  subcloud  down- 
draft  of  10  m/sec. 


z 

(m) 

(K? 

T 

(K) 

S 

a) 

M(xl0  3 ) 
(kg/m^) 

R  Z(xl05) 

(nm/hr)  ( mfi/m ?) 

R  - 

25  mm/ hr 

1500 

278.0 

278.0 

100.0 

1.24 

25.0 

0.33 

1000 

282.9 

282.8 

76.7 

1.23 

24.8 

0.33 

500 

287.8 

287.5 

60.5 

1.18 

24.1 

0.32 

0 

292.7 

292.1 

48.8 

1.12 

23.0 

0.31 

R  - 

50  mm/hr 

1500 

278.0 

278.0 

100.0 

2.24 

50.0 

0.92 

1000 

282.9 

282.8 

77.0 

2.22 

49.6 

0.91 

500 

287.8 

287.4 

61.5 

2.16 

48.6 

0.90 

0 

292.7 

291.8 

50.5 

2.07 

46.9 

0.88 

R  = 

75  ram/hr 

1500 

278.0 

278.0 

100.0 

3.28 

75.0 

1.77 

1000 

282.9 

282.8 

77.4 

3.25 

74.5 

1.76 

500 

287.8 

287.3 

62.4 

3.18 

73.1 

1.74 

0 

292.7 

291.6 

52.0 

3.07 

71.0 

1.71 

R  - 

100  mm/hr 

1500 

278.0 

278.0 

100.0 

4.03 

100.0 

2.51 

1000 

282.9 

282.7 

77.6 

3.99 

99.4 

2.50 

500 

287.8 

287.2 

63.0 

3.91 

97.8 

2.47 

0 

292.7 

291.4 

52.9 

3.79 

95.4 

2.43 
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Table  5.  Vertical  distributions  of  thermodynamic  and 

hydrometeoric  quantities  for  Case  A  in  a  subcloud 
downdraft  of  15  m/sec. 


z 

(m) 

To 

(K) 

T 

(K) 

S  M(xl03) 

(%)  (kg/m3) 

R  Z(xl0~5) 

(mm/hr) (mm^/m3) 

R  = 

25  mm/hr 

1500 

278.0 

278.0 

100.0 

1.24 

25.0 

0.33 

1000 

282.9 

282.9 

76.5 

1.23 

24.8 

0.33 

500 

287.8 

287.6 

59.8 

1.19 

24.2 

0.32 

0 

292.7 

292.3 

47.8 

1.14 

23.3 

0.31 

R  = 

50  mm/hr 

1500 

278.0 

278.0 

100.0 

2.24 

50.0 

0.92 

1000 

282.9 

282.8 

76.6 

2.22 

49.7 

0.92 

500 

287.8 

287.5 

60.5 

2.17 

48.7 

0.90 

0 

292.7 

292.1 

48.9 

2.09 

47.3 

0.88 

R  = 

75  mm/hr 

1500 

273.0 

278.0 

100.0 

3.28 

75.0 

1.77 

1000 

282.9 

282.8 

76.9 

3.25 

74.6 

1.76 

500 

287.8 

287.4 

61.1 

3.19 

73.3 

1.74 

0 

292.7 

291.9 

49.9 

3.09 

71.4 

1.71 

R  = 

100  mm/hr 

1500 

278.0 

278.0 

100.0 

4.03 

100.0 

2.51 

1000 

282.9 

282.8 

77.0 

4.00 

99.5 

2.50 

500 

287.8 

287.4 

61-5 

3.92 

98.1 

2.47 

0 

292.7 

291.8 

50.6 

3.82 

95.9 

2.43 

44 


for  constant  subcloud  downdrafts  of  5,  10,  and  15  m/sec,  respectively. 
The  subcloud  temperatures  resulting  from  a  dry-adiabatic  lapse  rate 
(Tp)  are  included  in  the  tables  for  comparison  with  the  computed  rain- 
shaft  temperature  (T) . 

The  details  of  the  vertical  variation  of  relative  humidity  are 
shown  in  Figs.  4  and  5  for  different  values  of  downdraft  and  for 
precipitation  rates  of  25  and  100  mm/hr,  respectively.  Clearly  the 
lowest  subcloud  relative  humidity  is  found  in  the  downdraft  with  the 
smallest  precipitation  rate  and  highest  downdraft  velocity.  The 
relative  humidity  throughout  the  subcloud  layer  should  approach  a 
value  of  100%  as  the  speed  of  the  downdraft  approaches  zero,  and  a 
dry-adiabatic  profile  as  the  speed  of  the  downdraft  becomes  very  large. 
The  relative  humidity  distributions  are  determined  by  the  evaporation 
of  precipitation,  which  is  affected  by  two  factors:  first,  the  smaller, 
slower-moving  drops  have  a  longer  time  available  to  evaporate  and  thus 
help  attain  a  higher  relative  humidity;  and,  second,  at  higher 
precipitation  rates  the  larger  concentration  of  small  drops  results 
in  more  evaporation  and  higher  relative  humidity. 

Changes  in  drop  radius  due  to  evaporation  are  relatively  large 
for  smaller  drops  (see  Appendix  C) .  One  way  to  see  this  more  closely 
is  through  the  evaporation  time  constant  which  is  the  time  required 
for  a  drop  to  be  reduced  by  evaporation  to  1/e  of  its  initial  mass. 

Time  constants  were  calculated  for  several  values  of  relative 
humidities  by  using  the  definition 


T 

r 


(70) 


HEIGHT  (m) 
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5. 


Same  as  in  Fig.  4  except  for  an  initial  precipitation 
rate  of  100  nnn/hr. 
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where  e  is  the  base  of  Naperian  logarithms  and  is  the  time  constant 
In  addition,  dMr/dt  is  computed  from  (C.2).  Fig.  6  illustrates 
the  results  of  the  time  constant  calculations.  As  the  environmental 
relative  humidity  decreases,  the  time  constant  for  the  small  rain¬ 
drops  decreases  to  the  point  where  very  few  of  these  drops  would 
survive  falling  through  one  height  interval  of  Az  =  25  m.  For 
this  reason,  drop-size  distributions  with  a  large  number  of  smaller 
drops  falling  at  a  slower  rate  give  relative  humidity  values  nearer  to 
saturation  than  with  distributions  dominated  by  larger,  faster-falling 
raindrops.  These  conclusions  are  consistent  with  the  results  of 
Kamburova  and  Ludlam  (1966),  who  concluded  that  raindrop  size  had  an 
appreciable  influence  on  evaporation,  which  proceeds  more  efficiently 
with  smaller  drop  sizes.  These  conclusions  also  are  consistent  with 
those  of  Das  and  Subbarao  (1972)  who  show  that  the  rainshaft  sub¬ 
saturation  is  maximum  with  smaller  water  content,  larger  drop  sizes, 
and  stronger  downdrafts. 

The  above  results  may  be  an  explanation  of  the  "humidity  dip"  ob¬ 
served  by  the  U.S.  Thunderstorm  Project  (Byers  and  Braham,  1949)  in 
strong  rain  shafts  accompanying  thunderstorms. 

Byers  and  Braham  suggested  that  one  reason  for  downdraft  subsaturation 
was  an  evaporation  rate  too  slow  to  provide  for  the  increase  in 
saturation  mixing  ratio  as  the  air  descended.  In  this  case,  the  air 
of  the  downdraft  would  be  heated  at  a  rate  between  the  moist  and  the 
dry-adiabatic  lapse  rates. 

Rainfall  rate  and  radar  reflectivity  are  not  significantly 
affected  by  evaporation  alone,  as  can  be  seen  from  Tables  3-5. 
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Figs.  7  and  8  depict  the  drop-size  distribution  1500  m  below 
cloud  base  for  the  conditions  presented  in  Figs.  4  and  5,  respectively. 
The  increased  loss  of  drops  at  the  small  radius  end  is  quite  obvious 
in  these  figures.  On  the  other  hand,  the  number  concentrations  of  drops 
larger  than  0.15  cm  in  radius  are  seen  to  change  very  little  thereby 
implying  that  the  number  of  drops  leaving  each  drop  category  due  to 
evaporation  is  nearly  equal  to  the  number  of  drops  entering  from  the 
next  larger  category. 

An  intuitively  apparent  result  is  seen  in  Figs.  7  and  8  in  which 
the  peaks  in  number  concentration  are  shifted  to  smaller  radii  for 
stronger  downdraft  velocities  and  higher  initial  precipitation  rates. 

It  is  easy  to  visualize  from  these  figures  that  in  the  limit 
as  the  downdraft  approaches  very  high  values,  the  peak  in  the  drop- 
size  distribution  will  continue  to  increase  and  shift  toward  small 
radii  until  it  mirrors  the  cloud-base  distribution.  The  data  used  to 
plot  Figs.  7  and  8  are  contained  in  Tables  6  and  7,  respectively.  The 
number  concentrations  of  drops  at  the  surface  are  maximum  for  large 
initial  rainfall  rates  and  strong  downdraft  velocities. 

The  drop-size  distributions  obtained  by  Kintigh  and  Das  (1970) 
for  pure  evaporation  showed  the  peak  near  a  radius  of  0.01  cm;  however, 
their  results  were  characterized  by  several  empty  drop-size  classes. 
Brazier-Smith  eit  al.  (1973)  found  that  the  number  concentrations  of 
all  drop  sizes  were  reduced  by  evaporation  with  a  preferential  re¬ 
duction  at  the  small-radius  end  of  the  spectrum.  Their  results  did 
not  show  a  peak  in  the  number  concentration  at  0.015-0.020  cm  radius 
because  they  assumed  a  constant  relative  humidity  gradient,  which  only 
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Radius  (cm) 


Drop-size  distributions  1500  m  below  cloud  base  for 
various  constant  downdrafts.  The  cloud  base  drop- 
size  distribution  is  Marshall-Palmer  for  R  =  25  mm/hr 
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Table  6.  Drop-size  distributions  at  cloud  base  (Marshall- 
Palmer  distribution  for  a  precipitation  rate  of 
25  mm/hr)  and  1500  m  below  cloud  base  for  constant 
downdrafts  of  5,  10,  and  15  m/sec  for  Case  A. 


-3  -1 

Drop  Concentration  (m  cm  ) 


Radius 

Cloud 

1500  m 

Below  Cloud 

Base 

(cm) 

Base 

5  m/sec 

10  m/sec 

15  m / sec 

0.0042 

1.340xl05 

1. 890x1 04 

3.015xl04 

3.999xl04 

0.0048 

1.311xl05 

2 . 189xl04 

3.527xl04 

4.648X104 

0.0053 

1.280x105 

2. 243xl04 

3.567xl04 

4.657xl04 

0.0060 

1.245xl05 

2.242xl04 

3.523xl04 

4.557xl04 

0.0067 

1.208x105 

2. 292xl04 

3.555xl04 

4.553xl04 

0.0076 

1. 167xl05 

2. 376xl04 

3 . 630xl04 

4.601xl04 

0.0085 

1 . 123xl05 

2.494xl04 

3.743x10^ 

4.689xl04 

0.0095 

1.075x105 

2. 644x1 04 

3.882xl04 

4.802xl04 

0.0107 

1.024xl05 

2.812xl04 

4 . 040xl04 

4.932xl04 

0.0120 

9.692xl04 

2. 998xl04 

4 . 202xl04 

5.062xl04 

0.0135 

9.115xl04 

3. 189xl04 

4.350xl04 

5 . 163xl04 

0.0151 

8.508xl04 

3.362X104 

4.445xl04 

5 . 209xl04 

0.0170 

7.875xl04 

3.400xl04 

4.383xl04 

5 . 048xl04 

0.0190 

7 . 220xl04 

3.479xl04 

4.331xl04 

4.901xl04 

0.0214 

6.550xl04 

3.463X104 

4.179X104 

4.657xl04 

0.0240 

5.871xl04 

3.394xl04 

3.974xl04 

4.362xl04 

0.0269 

5 . 193xl04 

3 . 246xl04 

3 . 698xl04 

4.003xl04 

0.0302 

4.525xl04 

3.038xl04 

3.379x10^ 

3.615X104 

0.0339 

3.876xl04 

2.757xl04 

2.999xl04 

3.175xl04 

0.0381 

3. 259xl04 

2.438xl04 

2.605xl04 

2. 733xl04 

0.0427 

2.682xl04 

2.099xl04 

2. 205x1 04 

2 . 294xl04 

0.0479 

2. 155xl04 

1 . 757xl04 

1.821xl04 

1.881xl04 

0.0538 

1.686xl04 

1.420xl04 

1.456xl04 

1.495xl04 

0.0604 

1 . 280xl04 

1. 106xl04 

1.126xl04 

1.150xl04 

0.0678 

9 . 392xl03 

8 . 311xl03 

8.407x10;* 

8.552xl03 

0.0761 

6.637x10;* 

5 . 991xl03 

6.033x10, 

6. 116xl03 

0.0854 

4.495xl03 

4 . 126xl03 

4. 140xl03 

4. 184xl03 

0.0959 

2. 902xl03 

2.700xl03 

2.703xl03 

2. 726xl03 

0.1076 

1.776xl03 

1.672x10, 

1.671xl03 

1.682xl03 

0.1208 

1.024x10, 

9 .721x10, 

9.701x10^ 

9 . 750x10^ 

0.1356 

5.513x10, 

5 . 263x10, 

5.247x10, 

5.269xl02 

0.1522 

2.753x10" 

2.643x10, 

2.634x10, 

2.642x10, 

0.1708 

1.263xl02 

1 . 215x10, 

1.210xl02 

1.213xl0f 

0.1918 

5.265X101 

5 .077x10, 

5.054x10* 

5 . 065x10, 

0.2152 

1.972x10*" 

1.904x10* 

1.895x10* 

1.899x10 

0.2416 

6 . 549x10® 

6.318x10° 

6 . 285x10° 

6.296x10° 

0.2712 

0.3044 

0.3417 

0.3835 

1. 900x10, 
4.739x10, 
9.970x10, 
1.733x10 

1.842x10, 

4.605x10, 

9.721x10, 

1.695x10^ 

1.833x10, 

4.584x10, 

9.682x10, 

1.689x10 

1.836x10, 

4.590x10, 

9.693x10, 

1.691x15 

i 
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Table  7, 


Radius 

(cm) 

0.0042 

0.0048 

0.0053 

0.0060 

0.0067 

0.0076 

0.0085 

0.0095 

0.0107 

0.0120 

0.0135 

0.0151 

0.0170 

0.0190 

0.0214 

0.0240 

0.0269 

0.0302 

0.0339 

0.0381 

0.0427 

0.0479 

0.0538 

0.0604 

0.0678 

0.0761 

0.0854 

0.0959 

0.1076 

0.1208 

0.1356 

0.1522 

0.1708 

0.1918 

0.2152 

0.2416 

0.2712 

0.3044 

0.3417 

0.3835 


Drop-size  distributions  at  cloud  base  (Marshall- 
Palmer  distribution  for  a  precipitation  rate  of 
100  mm/hr)  and  1500  m  below  cloud  base  for  constant 
downdrafts  of  5,  10,  and  15  m/sec  for  Case  A. 


Cloud 

Base 


Drop  Concentration  (m"3cm“l) 

1500  m  Below  Cloud  Base 
5  m/sec  10  m/sec  15  m/sec 


4.826xl04 
5.617x10^ 
5.653x107 
5 . 540x10? 
5.544x10? 
5.617x10? 
5.743x10? 
5.907x10? 
6 . 103x10? 
6 . 309x10? 
6.494x10? 
6.626x10? 
6.506x10? 
6 .416x10, 
6 . 211x10? 
5 . 948x10? 
5.603x10? 
5.215x10? 
4 . 744x10? 
4.251x10? 
3.736x10? 
3.228x10, 
2.723x10? 
2.240x10? 
1.798x10, 
1.401x10? 
1.057x10 
7.682x10 
5.363x10 
3.574x10 
2.259x10 
1.351x10 
7.563x10 
3.945x10 
1.899x10 
8.330x10 
3.326x10 
1.184x10 
3.716x10 
1.012x10 


1.402xl05 

2.607xl04 

3.813x10^ 

1.380x105 

3.074x10? 

4 .465x10? 

1.355x10^ 

3.167x10? 

4.537x10? 

1.327x10? 

3.159x10? 

4.483x10? 

1.297x10^ 

3.225x10? 

4.526x10^ 

1 . 265x10]? 

3.340x10? 

4.629x10? 

1.229x10? 

3.504x10? 

4.782x10? 

1.190x10^ 

3.715x10? 

4.976x10? 

1.147x10^ 

3.956x10? 

5.199x10? 

1.101x10^ 

4.226x10^ 

5.439x10^ 

1.052x10^ 

4.508x10? 

5.670x10? 

9 . 994x10? 

4.770x10? 

5.859x10? 

9.435x10? 

4.845x10? 

5.824x10? 

8.844x10? 

4.989x10? 

5.829x10? 

8.225x10? 

5.010x10? 

5.715x10? 

7.581x10? 

4.971x10? 

5 . 540x10, 

6.918x10? 

4.835x10? 

5.276x10? 

6.243x10? 

4.625x10? 

4.958x10? 

5.564x10? 

4.311x10? 

4.547x10? 

4.889x10? 

3.938x10? 

4.103x10? 

4.228x10? 

3.526x10? 

3.630x10? 

3.592x10? 

3.095x10? 

3 . 154x10? 

2.992x10? 

2.640x10? 

2.673x10? 

2.436x10? 

2.191x10? 

•  2.207x10? 

1.935x10? 

1.770x10? 

1 . 776x10? 

1.494x10? 

1.387x10? 

1.388x10, 

1.117x10? 

1.050x10? 

1.049x10* 

8.065x10? 

7.657x10? 

7.639x10? 

5.594x10? 

5 . 358x10:: 

5 . 341x10? 

3 . 710x10? 

3.577x10? 

3.563x10^ 

2. 340x10? 

2.263x10? 

2.254x10? 

1.395x10? 

1 . 355x10? 

1 . 349x10? 

7.803x10? 

7.589x10? 

7.552x10, 

4.066x10? 

3.960x10? 

3.941x10, 

1.956x10? 

1.907x10^ 

1.897x10^ 

8.604x10? 

8.376x10? 

8.325x10? 

3.422x10? 

3.343x10? 

3.324x10? 

1.216x10,: 

1.190x10? 

1.183x10.: 

3.806x10° 

3.732x10° 

3.715x10° 

1.033x10 

1.016x10 

1.012x10° 
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allowed  a  70%  relative  humidity  1500  m  below  cloud  base.  As  Indicated 
in  Fig.  3,  the  subsaturation  in  a  downdraft  attains  much  greater  magni¬ 
tudes,  thereby  substantially  reducing  the  evaporation  time  constant  for 
the  small  drops. 

An  unexpected  feature  of  the  present  computations  is  an 
oscillation  in  the  drop  concentration  between  the  radii  of  0.005  and 
0.009  cm,  that  can  be  seen  in  Tables  6  and  7.  An  investigation  (see 
Appendix  D)  into  this  anomaly  showed  its  origin  to  be  in  the  evapora¬ 
tion  data  of  Kinzer  and  Gunn  (1951).  At  any  rate,  this  anomaly  has 
little  significance  to  the  main  thrust  of  the  present  study. 

b.  Case  B:  Evaporation  and  Collision-Coalescence 

Collision-coalescence  was  added  to  Case  A  using  (56)  and  (61). 

The  combined  effects  of  evaporation  and  collision-coalescence  are 
presented  in  Tables  8,  9,  and  10,  which  show  many  interesting  changes 
from  Case  A.  On  comparing  Tables  8-10  with  Tables  3-5,  the  following 
differences  become  apparent: 

1.  The  precipitation  rates  for  Case  B  increase  with  depth  below 
the  cloud  base  because  as  the  larger  drops  sweep  out  the  smaller  ones, 
the  liquid  water  represented  by  the  sum  of  their  masses  has  increased 
their  average  terminal  velocity  and  their  contribution  to  the  precipi¬ 
tation  rate.  This  is  generally  true  except  for  small  deviations  in 
the  case  of  25-  and  50-mm/hr  (initial)  precipitation  rates.  The  number 
of  larger  drops  initially  available  for  these  smaller  precipitation 
rates  is  much  less  than  for  precipitation  rates  of  75  and  100  mm/hr, 
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Table  8.  Vertical  distributions  of  thermodynamic  and 

hydrometeoric  quantities  for  Case  B  in  a  sub¬ 
cloud  downdraft  of  5  m/sec. 


z 

(m) 

(K) 

T 

(K) 

S 

(%) 

M(xl03) 

(kg/m3) 

R  Z(xl0"5) 

(mm/hr)  (mm^/tn^) 

R  *= 

25  mm/hr 

1500 

278.0 

278.0 

100.0 

1.24 

25.0 

0.33 

1000 

282.9 

282.8 

77.3 

1.22 

25.6 

0.36 

500 

287.8 

287.3 

61.9 

1.19 

25.7 

0.38 

0 

292.7 

291.8 

50.9 

1.14 

25.4 

0.40 

R  * 

50 

1500 

278.0 

278.0 

100.0 

2.24 

50.0 

0.92 

1000 

282.9 

282.7 

77.8 

2.22 

52.1 

1.05 

500 

287.8 

287.2 

63.3 

2.17 

53.2 

1.16 

0 

292.7 

291.4 

53.0 

2.12 

53.5 

1.26 

R  * 

75  mm/hr 

1500 

287.0 

287.0 

100.0 

3.28 

75.0 

1.77 

1000 

282.9 

282.7 

78.3 

3.25 

79.1 

2.09 

500 

287.8 

287.0 

64.4 

3.21 

81.5 

2.38 

0 

292.7 

291.2 

54.7 

3.15 

82.7 

2.64 

R  * 

100  mm/hr 

1500 

287.0 

287.0 

100.0 

4.03 

100.0 

2.51 

1000 

282.9 

282.6 

78.6 

4.00 

105.9 

3.02 

500 

287.8 

286.9 

65.1 

3.95 

109.2 

3.49 

0 

292.7 

291.0 

55.7 

3.89 

111.1 

3.91 
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Table  9.  Vertical  distributions  of  thermodynamic  and 

hydrometeoric  quantities  for  Case  B  in  a  sub¬ 
cloud  downdraft  of  10  m/sec. 


z 

(m) 

Td 

<K) 

T 

(K) 

S 

<%) 

M(xl03) 

(kg/m3) 

R  Z(xl0*5) 

(mm/hr)  (mm^ /m3  ) 

R  = 

25  mm/hr 

1500 

278.0 

278.0 

100.0 

1.24 

25.0 

0.33 

1000 

282.9 

282.8 

76.7 

1.23 

25.3 

0.35 

500 

287.8 

287.5 

60.4 

1.19 

25.2 

0.36 

0 

292.7 

292.2 

48.5 

1.14 

24.7 

0.37 

R  = 

50  mm/hr 

1500 

278.0 

278.0 

100.0 

2.24 

50.0 

0.92 

1000 

282.9 

282.8 

77.0 

2.22 

51.2 

1.00 

500 

287.8 

287.4 

61.2 

2.17 

51.7 

1.07 

0 

292.7 

292.0 

49.7 

2.11 

51.5 

1.13 

R  = 

75  mm/hr 

1500 

278.0 

278.0 

100.0 

3.28 

75.0 

1.77 

1000 

282.9 

282.8 

77.2 

3.25 

77.6 

1.98 

500 

287.8 

287.4 

61.8 

3.22 

79.0 

2.17 

0 

292.7 

291.8 

50.7 

3.13 

79.4 

2.33 

R  - 

100  mm/hr 

1500 

278.0 

278.0 

100.0 

4.03 

100.0 

2.51 

1000 

282.9 

282.8 

77.4 

4.00 

103.8 

2.85 

500 

287.8 

287.3 

62.2 

3.94 

105.9 

3.16 

0 

292.7 

291.7 

51.4 

3.86 

106.9 

3.44 
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Tahle  10.  Vertical  distributions  of  thermodynamic  and 

hydrometeoric  quantities  for  Case  B  in  a  sub¬ 
cloud  downdraft  of  15  m/sec. 


z 

(m) 

TD 

(K) 

T 

(K) 

CO  ^5 

M(xl03) 

(kg/ra3) 

R  Z(xl0-5) 

(mm/hr)  (mn^/m3) 

R  = 

25  mm/hr 

1500 

278.0 

278.0 

100.0 

1.24 

25.0 

0.33 

1000 

282.9 

282.9 

76.4 

1.23 

25.2 

0.34 

500 

287.8 

287.6 

59.8 

1.20 

25.0 

0.35 

0 

292.7 

292.3 

47.6 

1.15 

24.5 

0.36 

R  = 

50  mm/hr 

1500 

278.0 

278.0 

100.0 

2.24 

50.0 

0.92 

1000 

282.9 

282.8 

76.7 

2.22 

50.9 

0.98 

500 

287.8 

287.5 

60.3 

2.18 

51.1 

1.03 

0 

292.7 

292.2 

48.5 

2.12 

50.8 

1.07 

R  = 

75  mm/hr 

1500 

278.0 

278.0 

100.0 

3.28 

75.0 

1.77 

1000 

282.9 

282.8 

76.8 

3.26 

76.9 

1.92 

500 

287.8 

287.5 

60.8 

3.20 

77.8 

2.06 

0 

292.7 

292.0 

49.2 

3.13 

77.9 

2.18 

R  = 

100  mm/hr 

1500 

278.0 

278.0 

100.0 

4.03 

100.0 

2.51 

1000 

282.9 

282.8 

77.0 

4.00 

102.7 

2.76 

500 

287.8 

287.4 

61.1 

3.94 

104.3 

2.99 

0 

292.7 

292.0 

49.7 

3.86 

104.9 

3.20 

so  the  shift  of  liquid  water  from  the  smaller  to  the  larger  drop 
radii  is  not  as  great.  As  a  result,  the  precipitation  rate  increases 
initially,  but  then  decreases  as  the  evaporation  begins  to  dominate 
the  collision-coalescence  effect.  7  «  effect  of  evaporation  is  evident 
in  Tables  8-10,  because  total  liquid  water  decreases  from  cloud  base  to 
ground . 

2.  Radar  reflectivity,  like  precipitation  rate,  increases  with 
depth  below  cloud  base.  The  larger  drops  backscatter  a  higher  amount 
of  radar  energy,  so  even  though  the  total  amount  of  liquid  water  de¬ 
creases,  the  shift  of  liquid  water  from  the  smaller  to  the  larger 
drops  causes  the  radar  reflectivity  to  increase. 

3.  The  surface  temperature  is  warmer  and  the  relative  humidity 
lower  for  Case  B  because  the  larger  drops  sweep  out  the  smaller  ones, 
the  more  efficient  evaporators. 

Relative  humidity  profiles  for  Cases  A  and  B  are  presented  in 
Figs.  9  and  10  for  initial  precipitation  rates  of  25  and  100  mm/hr, 
respectively.  Fig.  9  illustrates  that  little  difference  exists  be¬ 
tween  Cases  A  and  B  for  the  lower  precipitation  rate  and  strong  down- 
draft  velocity.  In  each  case,  the  profile  approaches  the  dry  adia¬ 
batic  curve  because  the  strong  downdraft  velocity  does  not  allow  enough 
time  for  evaporation.  Fig.  10  shows  that  as  the  number  of  larger 
drops  increases,  the  subsaturation  increases.  Thus  the  difference  of 
relative  humidity  between  Cases  A  and  B  is  maximum  for  large  initial 
precipitation  rate  and  weak  downdraft  velocity. 

The  surface  drop-size  distributions  resulting  from  Case  B  are 
presented  in  Figs.  11  and  12  for  initial  precipitation  rates  of  25  and 
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Fig.  10.  Same  as  Fig.  9  except  for  a  constant  downdraft  of  5  m/sec 
and  initial  precipitation  rate  of  100  nnn/hr. 
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100  mm/hr,  respectively.  The  effect  of  collision-coalescence  in 
depleting  the  number  of  small  drops  is  readily  apparent.  Fig.  12 
illustrates  the  sharp  drop  in  the  number  concentration  of  small  drops 
as  the  downdraft  velocity  decreases.  This  is  due  to  the  longer  time 
the  large  drops  have  to  sweep  out  the  smaller  ones. 

Figs.  13  and  14  present  comparisons  of  the  drop-size  distributions 
at  the  surface,  obtained  in  Cases  A  and  B.  In  Fig.  13  which  relates 
to  the  cloud-base  precipitation  rate  of  25  mm/hr,  the  number  of  small 
drops  decreases  significantly  from  Case  A  to  Case  B,  but  this  decrease 
is  even  sharper  in  Fig.  14  where  the  initial  concentration  of  larger 
drops  is  higher. 

c.  Case  C:  Evaporation,  Collision-Coalescence,  and  Aerodynamic 
Breakup 

The  computations  in  Case  C  were  performed  by  adding  (62)  and  (63) 
to  Case  B.  Vertical  distributions  of  thermodynamic  and  hydrometeor ic 
quantities  resulting  from  these  computations  are  given  in  Tables  11, 

12,  and  13. 

A  comparison  of  Tables  11-13  with  Tables  8-10  points  up  one  fact 
quite  clearly:  there  is  little  difference  between  the  two  cases 
for  the  initial  precipitation  rate  of  25  mm/hr,  but  the  results  show 
marked  difference  as  the  initial  precipitation  rate  increases.  Rela¬ 
tive  humidity  for  25  mm/hr  is  lower  at  the  surface  for  all  downdraft 
velocities  in  Case  C  when  compared  to  Case  B,  but  by  only  0.1  percent 
or  less.  This  difference  increases  to  about  0.3  percent  for  the 
intensity  of  100  mm/hr.  The  most  significant  difference  in  Cases  B 


Radius  (cm) 


Drop-size  distributions  1500  m  below  cloud  base  for 
Case  A  and  Case  B,  a  constant  downdraft  of  5  m/sec 
and  an  initial  precipitation  rate  of  25  mm/hr. 
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Table  11.  Vertical  distributions  of  thermodynamic  and 

hydrometeoric  quantities  for  Case  C  in  a  sub¬ 
cloud  downdraft  of  5  m/sec. 


z 

(m) 

Td 

(K) 

T 

(K) 

S 

(7.) 

M(xl03) 

(kg/m3) 

R  Z(xl0"5) 

(mm/hr) (mm^/m^) 

R  = 

25  mm/hr 

1500 

278.0 

278.0 

100.0 

1.24 

25.0 

0.33 

1000 

282.9 

282.8 

77.3 

1.22 

25.7 

0.37 

500 

287.8 

287.3 

61.9 

1.19 

25.9 

0.40 

0 

292.7 

291.8 

50.8 

1.14 

25.6 

0.43 

R  = 

50 

1500 

278.0 

278.0 

100.0 

2.24 

50.0 

0.92 

1000 

282.9 

282.7 

77.8 

2.22 

52.9 

1.18 

500 

287.8 

287.2 

63.2 

2.18 

54.4 

1.39 

0 

292.7 

291.5 

52.7 

2.12 

55.0 

1.54 

R  = 

75  mm/hr 

1500 

278.0 

278.0 

100.0 

3.28 

75.0 

1.77 

1000 

282.9 

282.7 

78.2 

3.25 

81.5 

2.58 

500 

287.8 

287.1 

64.1 

3.21 

84.4 

3.09 

0 

292.7 

291.2 

54.3 

3.15 

85.0 

3.36 

R  = 

100  mm/hr 

1500 

278.0 

278.0 

100.0 

4.03 

100.0 

2.51 

1000 

282.9 

282.7 

78.5 

4.00 

109.4 

3.91 

500 

287.8 

287.0 

64.7 

3.95 

112.7 

4.58 

0 

292.7 

291.1 

55.4 

3.89 

112.9 

4.83 

V 
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Table  12.  Vertical  distributions  of  thermodynamic  and 

hydrometeoric  quantities  for  Case  C  in  a  sub¬ 
cloud  downdraft  of  10  m/sec. 


z 

(in) 

Td 

(K) 

T 

(K) 

S 

<Z) 

MCxlO3) 

(kg/m3) 

R  Z(xl0-5) 

(mm/hr) (mm^/m3) 

R  = 

25  mm/hr 

1500 

278.0 

278.0 

100.0 

1.24 

25.0 

0.33 

1000 

282.9 

282.8 

76.7 

1.23 

25.4 

0.36 

500 

287.8 

287.5 

60.3 

1.19 

25.3 

0.38 

0 

292.7 

292.2 

48.4 

1.14 

24.9 

0.39 

R  = 

50  mm/hr 

1500 

278.0 

278.0 

100.0 

2.24 

50.0 

0.92 

1000 

282.9 

282.8 

77.0 

2.22 

51.9 

1.10 

500 

287.8 

287.4 

61.1 

2.17 

52.8 

1.26 

0 

292.7 

292.0 

49.6 

2.11 

53.0 

1.37 

R  = 

75  mm/hr 

1500 

278.0 

278.0 

100.0 

3.28 

75.0 

1.77 

1000 

282.9 

282.8 

77.2 

3.25 

79.8 

2.37 

500 

287.8 

287.4 

61.6 

3.20 

82.0 

2.81 

0 

292.7 

291.8 

50.4 

3.14 

82.6 

3.07 

R  = 

100  mm/hr 

1500 

278.0 

278.0 

100.0 

4.03 

100.0 

2.51 

1000 

282.9 

282.8 

77.3 

4.00 

107.0 

3.58 

500 

287.8 

287.3 

62.0 

3.94 

110.2 

4.23 

0 

292.7 

291.8 

51.0 

3.87 

110.6 

4.54 
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Table  13.  Vertical  distributions  of  thermodynamic  and 

hydrometeoric  quantities  for  Case  C  in  a  sub¬ 
cloud  downdraft  of  15  m/sec. 


z 

(m) 

td 

(K) 

T 

(K) 

S 

(%) 

MCxlO3) 

(kg/m3) 

R  Z(xl0-5) 

(mm/hr)  (mnrVm3) 

R  = 

25  mm/hr 

1500 

278.0 

278.0 

100.0 

1.24 

25.0 

0.33 

1000 

282.9 

282.9 

76.4 

1.23 

25.3 

0.35 

500 

287.8 

287.6 

59.8 

1.20 

25.1 

0.36 

0 

292.7 

292.3 

47.5 

1.15 

24.7 

0.37 

R  = 

50  mm/hr 

1500 

278.0 

278.0 

100.0 

2.24 

50.0 

0.92 

1000 

282.9 

282.8 

76.7 

2.22 

51.4 

1.06 

500 

287.8 

287.1 

61.7 

2.18 

52.1 

1.18 

0 

292.7 

292.2 

48.4 

2.12 

52.1 

1.28 

R  = 

75  mm/hr 

1500 

278.0 

278.0 

100.0 

3.28 

75.0 

1.77 

1000 

282.9 

282.8 

76.8 

3.26 

78.5 

2.25 

500 

287.8 

287.5 

60.7 

3.21 

80.5 

2.62 

0 

292.7 

292.1 

49.0 

3.14 

82.2 

2.87 

R  = 

100  mm/hr 

1500 

278.0 

278.0 

100.0 

4.03 

100.0 

2.51 

1000 

282.9 

282.8 

76.9 

4.00 

105.6 

3.38 

500 

287.8 

287.5 

61.0 

3.95 

108.5 

3.97 

0 

292.7 

292.0 

49.4 

3.88 

109,2 

4.31 

and  C  subcloud  quantities  is  the  radar  reflectivity  and  rainfall  rate. 
Both  of  these  quantities  increase  due  to  the  aerodynamic  breakup  of  the 
largest  drops  into  smaller  drop  categories. 

In  order  to  have  a  better  understanding  of  the  dependence  of  aero¬ 
dynamic  breakup  on  drop  radius,  (28)  and  (29)  were  investigated  in 
more  detail.  Table  14  shows  how  P^,  the  probability  of  drop  breakup 
per  unit  time,  varies  with  drop  radius.  It  is  seen  that  only  the 

Table  14.  Probability  of  aerodynamic 
drop  breakup. 


Radius 

(cm) 

Pk  _ 
(sec 

0.2416 

0.001 

0.2712 

0.003 

0.3044 

0.009 

0.3417 

0.033 

0.3835 

0.135 

largest  category  of  drops  has  a  significant  probability  of  aerodynamic 
breakup.  Aerodynamic  breakup,  consequently,  has  a  greater  effect  on 
drop-size  distributions  having  large  initial  concentrations  of  large 
drops  which  will  be  the  case  for  Marshall-Palmer  distributions  at  a  high 
precipitation  intensity.  The  fragment  probability  given  by  Q(k,  k’)  in 
(29)  determines  how  many  fragment  drops  fall  in  each  smaller  radius 
category  k  due  to  breakup  of  a  large  drop  of  radius  k',  and  shows  that 
the  fragment  population  will  be  maximum  at  the  drop  radius  r^/7- 
Since  the  mid-point  radius  in  the  largest  category  is  r^,  =  0.3835  cm, 
the  fragment  maximum  will  be  at  a  radius  of  about  0.055  cm.  This 
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fragment  maximum  will  decrease  with  decreasing  r^, ,  so  when  a  drop 
distribution  has  a  high  concentration  of  large  drops,  the  small  radius 
end  of  the  distribution  will  benefit  most  from  aerodynamic  breakup. 

Figs.  15  and  16  show  the  drop-size  distributions  resulting  from 
Case  C,  while  a  comparison  of  drop-size  distributions  for  Cases  B  and 
C  is  shown  in  Figs.  17  and  18.  The  results  in  Fig.  17  show  little 
difference  between  Cases  B  and  C;  however,  the  contribution  to  the 
concentration  of  small  drops  due  to  breakup  of  the  largest  drops  is 
apparent.  The  data  used  to  plot  Figs.  17  and  18  are  given  in  Tables 
15  and  16,  respectively.  Note,  in  Table  16,  that  the  drop  concentra¬ 
tion  for  the  two  largest  radius  categories  is  smaller  for  Case  C  than 
B.  This  can  be  expected  from  the  probabilities  in  Table  14.  Fig.  18 
illustrates  the  dramatic  effect  of  a  high  concentration  of  large  drops. 
These  results  are  consistent  with  the  results  of  Srivastava  (1971),  who 
noted  that  aerodynamic  breakup  has  a  tendency  to  flatten  the  equilib¬ 
rium  distribution  compared  to  the  distributions  observed  in  nature. 

Aerodynamic  breakup  has  little  effect  on  drop-size  distributions 
having  few  large  drops.  For  that  reason,  this  microphysical  process 
is  likely  to  be  important  in  cold  rather  than  warm  rain  processes. 
Srivastava  (1971)  suggested  that  collisional  breakup  may  help  produce 
better  agreement  between  theory  and  observation. 

d.  Case  D:  Evaporation,  Collision-Coalescence,  and  Collisional 

Breakup 

Case  D  includes  the  microphysical  processes  of  evaporation, 
collision-coalescence  with  the  coalescence  efficiency  of  Whelpdale 
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Same  as  Fig.  7  except  for  Case  C 
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Fig.  16.  Same  as  Fig.  8  except  for  Case  C. 
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Cloud  Base,  R  =  25  mm/hr 
Case  B 


Radius  (cm) 


Fig.  17.  Drop-size  distributions  1500  m  below  cloud  base  for 
Case  B  and  Case  C,  a  constant  downdraft  of  5  m/sec 
and  an  initial  precipitation  rate  of  25  mm/hr. 
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Table  15.  Drop-size  distributions  at  cloud  base  (Marshall- 

Palmer)  and  1500  m  below  cloud  base  for  Case  B  and 
Case  C.  The  constant  downdraft  is  5  m/sec  and  the 
initial  precipitation  rate  is  25  ram/hr. 


-3 

-1 

Drop 

Concentration  (m 

cm  ) 

Radius 

Cloud 

1500  m  Below 

Cloud  Base 

(cm) 

Base 

Case  B 

Case  C 

0.0042 

1.340x10 

5 

q 

5.739x10  I 

5.726x10 

3 

0.0048 

1.311x10 

q 

5.493x10  , 

7.100x10 

J 

0.0053 

1.280x10 

q 

5.658x10  i 

7.222x10 

3 

0.0060 

1.245x10 

5.697x10  " 

7.178x10 

J 

0.0067 

1.208x10 

j 

q 

5.873x10  , 

7.299x10 

J 

0.0076 

1.167x10 

q 

6.148x10  - 

7.528x10 

J 

0.0085 

1.123x10 

q 

6.526x10  ^ 

7.868x10 

J 

0.0095 

1.075x10 

q 

7.011x10  , 

8.317x10 

3 

0.0107 

1.024x10 

J 

1 V 

7.577x10  ^ 

8.840x10 

3 

0.0120 

9.692x10 

A 

8.233x10  " 

9.444x10 

J 

0.0135 

9.115x10 

A 

8.969x10  ^ 

1.012x10 

4 

A 

0.0151 

8.508x10 

A 

9.731x10  "T 

1.079x10 

0.0170 

7.875x10 

A 

1.020x10  7 

1.114x10 

4 

0.0190 

7.220x10 

A 

1.092x10  J 

1.173x10 

4 

0.0214 

6.550x10 

4 

A 

1.143x10  7 

1.210x10 

4 

0.0240 

5.871x10 

4 

A 

1.190x10 

1.242x10 

4 

0.0269 

5.193x10 

4 

A 

1.221x10  7 

1.258x10 

4 

A 

0.0302 

4.525x10 

4 

A 

1.239x10  7 

1.261x10 

4 

A 

0.0339 

3.876x10 

4 

1.235x10  7 

1.243x10 

0.0381 

3.259x10 

4 

A 

1.209x10  7 

1.205x10 

4 

0.0427 

2.682x10 

A 

1.170x10  7 

1.156x10 

4 

0.0479 

2.155x10 

4 

A 

1.103x10  7 

1.081x10 

4 

0.0538 

1.686x10 

A 

1.001x10  ^ 

9.750x10 

3 

0.0604 

1.280x10 

4 

8.714x10  ^ 

8.455x10 

J 

0.0678 

9.392x10 

J 

7.249x10  , 

7.015x10 

3 

0.0761 

6.637x10 

J 

5.736x10  J 

5.546x10 

J 

0.0854 

4.495x10 

j 

4.311x10  ' 

4.171x10 

3 

0.0959 

2.902x10 

J 

3.038x10  t 

2.949x10 

3 

0.1076 

1.776x10 

J 

1.998x10  -1 

1.951x10 

3 

0.1208 

1.024x10 

J 

o 

1.224x10  i 

1.207x10 

3 

0.1356 

5.513x10 

9 

6.893x10  , 

6.894x10 

Z 

0.1522 

2.753x10 

9 

3.544x10  , 

3.622x10 

Z 

0.1708 

1.263x10 

1.654x10 

1.745x10 

z 

1 

0.1918 

5.265x10 

1 

i 

6.930x10  | 

7.661x10 

0.2152 

1.972x10 

A 

2.594x10  * 

3.059x10 

i 

0.2416 

6.549x10 

U 

A 

8.531x10  ^ 

1.087x10 

1 

0.2712 

1.900x10 

2.465x10  j 

3.260x10 

u 

0.3044 

4.739x10' 

1  I 

,  9 

6.103x10": 

6.665x10" 

1 

0.3417 

9.970x10 

.  9 

1.277x10  : 

6.837x10' 

z 

0.3835 

1.733x10 

2.208x10' 

5.384x10' 

3 

\ 
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Table  16.  Drop-size  distributions  at  cloud  base  (Marshall- 

Palmer)  and  1500  m  below  cloud  base  for  Case  B  and 
Case  C.  The  constant  downdraft  is  5  m/sec  and  the 
initial  precipitation  rate  is  100  mm/hr. 


-3  -1 

Drop  Concentration  (m  cm  ) 


Radius 

(cm) 

Cloud 

Base 

1500  m  Below 

Case  B 

Cloud  Ba; 

Case  C 

0.0042 

1.402x10  5 

3.634x10  j 

8.313x10 

0 . 0048 

1.380x10  Z 

4.340x10  j 

9.507x10 

0.0053 

1.355x10  * 

4.539x10  j 

9.521x10 

0.0060 

1.327x10 

4.627x10  2 

9.246x10 

0.0067 

1.297x10  Z. 

4.840x10  2 

9.139x10 

0.0076 

1.265x10  Z 

5.155x10  j 

9.120x10 

0.0085 

1.229x10  Z 

5.586x10  f 

9.179x10 

0.0095 

1.190x10  Z 

6.150x10  2 

9.300x10 

0.0107 

1.147x10  Z 

6.844x10  j 

9.431x10 

0.0120 

1.101x10  z. 

7.704x10  2 

9.582x10 

0.0135 

1.052x10  , 

8.756x10  2 

9.741x10 

0.0151 

9.994x10  7 

9.995x10 

9.869x10 

0.0170 

9.435x10  7 

1.114x10  3 

9.704x10 

0  0190 

8.844x10  7 

1.283x10  3 

9.723x10 

0.0214 

8.225x10  7 

1.460x10  3 

9.577x10 

0.0240 

7.581x10  7 

1.676x10  3 

9.409x10 

0.0269 

6.918x10  , 

1.924x10  3 

9.173x10 

0.0302 

6.243x10  , 

2.223x10  3 

8.902x10 

0.0339 

5.564x10  , 

2.578x10  " 

8.600x10 

0.0381 

4.889x10  , 

2.989x10 

8.251x10 

0.0427 

4.228x10  7 

3.520x10 

7.931x10 

0.0479 

3.592x10  , 

4.087x10  3 

7.571x10 

0.0538 

2.992x10  , 

4.628x10  3 

7.151x10 

0.0604 

2.436x10 

5.094x10 

6.673x10 

0.0678 

1.935x10 

5.409x10  3 

6.125x10 

0.0761 

1.494x10  7 

5.550x10  3 

5.531x10 

0.0854 

1.117x10  Z. 

5.548x10  3 

4.928x10 

0.0959 

8.065x10 

5.161x10  3 

4.237x10 

0.1076 

5.594x10  Z. 

4.575x10  3 

3.553x10 

0.1208 

3.710x10 

3.800x10  3 

2.873x10 

0.1356 

2.340x10  3 

2.875x10  3 

2.201x10 

0.1522 

1.395x10  j 

1.977x10  3 

1.599x10 

0.1708 

7.803x10  j 

1.234x10  2 

1.094x10 

0.1918 

4.066x10  j 

6.865x10  2 

7.007x10 

0.2152 

1.956x10 

3.449x10  j 

4.210x10 

0.2416 

8.604x10  : 

1.518x10  \ 

2.294x10 

0.2712 

3.422x10 

6.040x10  ^ 

1.104x10 

0.3044 

1.216x10  „ 

2.109x10 

3.959x10 

0.3417 

3.806x10  : 

6.514x10  ^ 

7.896x10 

0.3835 

1.034x10 

1.746x10 

8.490x10 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

2 

2 

2 

2 

1 

0 

-1 
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and  List  (1971),  and  collisional  breakup  incorporating  the  experimental 
fragment  size  distributions  of  McTaggart-Cowan  and  List  (1975)  as 
formulated  by  List  and  Gillespie  (1976). 

Vertical  distributions  of  thermodynamic  and  hydrometeoric  quanti¬ 
ties  resulting  from  Case  D  for  initial  precipitation  rates  of  25,  50, 

75,  and  100  mm/hr  are  given  in  Tables  17,  18,  and  19  for  constant 
downdrafts  of  5,  10,  and  15  m/sec,  respectively. 

On  comparing  the  results  in  Case  D  as  detailed  in  Tables  17-19 
with  those  in  Case  B  (Tables  8-10)  and  Case  C  (Tables  11-13),  the 
following  points  become  clear. 

1.  The  difference  in  the  rainshaft  temperature  and  relative 
humidity  between  Cases  B  and  D  is  greater  for  a  weak  downdraft 

(5  m/sec)  than  for  stronger  downdrafts  (10  and  15  m/sec).  Collisional 
breakup  is  a  more  efficient  process  than  aerodynamic  breakup  in  deplet¬ 
ing  the  large  drops  and  creating  small  fragments.  With  weaker  down- 
drafts  the  small  drops  have  more  time  to  evaporate,  thus  yielding  a 
higher  rainshaft  relative  humidity  and  lower  rainshaft  temperature.  With 
stronger  downdrafts  the  small  fragments  are  not  given  as  much  time  to 
evaporate,  so  the  difference  between  Cases  B  and  D  is  less  pronounced. 

2.  Both  the  rainfall  rate  and  the  radar  reflectivity  are  greatly 
reduced  for  Case  D  when  compared  to  Cases  B  and  C.  This  is  due  to  the 
significant  loss,  in  Case  D,  of  large  drops  which,  with  their  high 
terminal  velocities,  contribute  heavily  to  the  overall  rainfall  rate. 
These  large  drops  also  are  efficient  scatterers  of  the  radar  signal, 
so  their  loss  results  in  a  significant  loss  of  radar  ref lectivity . 
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Table  17.  Vertical  distributions  of  thermodynamic  and 

hydrometeoric  quantities  for  Case  D  in  a  sub¬ 
cloud  downdraft  of  5  m/sec. 


z 

(m) 

td 

(K) 

T 

(K) 

S  M(xl03) 
(%)  (kg/m3) 

5 

R  Z(xl0-J) 

(mm/hr)  (tran^/cm3) 

R  = 

25  mm/hr 

1500 

278.0 

278.0 

100.0 

1.24 

25.0 

0.33 

1000 

282.9 

282.8 

77.4 

1.22 

24.1 

0.22 

500 

287.8 

287.3 

62.4 

1.18 

23.1 

0.19 

0 

292.7 

291.6 

51.7 

1.11 

22.0 

0.17 

R  = 

50  mm/hr 

1500 

278.0 

278.0 

100.0 

2.24 

50.0 

0.92 

1000 

282.9 

282.7 

78.1 

2.21 

46.8 

0.45 

500 

287.8 

287.0 

64.4 

2.15 

45.2 

0.38 

0 

292.7 

291.1 

55.1 

2.07 

43.7 

0.35 

R  = 

75  mm/hr 

1500 

278.0 

278.0 

100.0 

3.28 

75.0 

1.77 

1000 

282.9 

282.6 

78.8 

3.25 

71.3 

0.69 

500 

287.8 

286.8 

66.2 

3.17 

69.3 

0.58 

0 

292.7 

290.7 

58.1 

3.07 

67.4 

0.55 

R  = 

100  mm/hr 

1500 

278.0 

278.0 

100.0 

4.03 

100.0 

2.51 

1000 

282.9 

282.6 

79.2 

3.99 

90.3 

0.85 

500 

287.8 

286.6 

67.4 

3.91 

88.0 

0.73 

0 

292.7 

290.4 

60.1 

86.0 

86.0 

0.69 
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Table  18.  Vertical  distribution  of  thermodynamic  and 

hydrometeorlc  quantities  for  Case  D  in  a  sub¬ 
cloud  downdraft  of  10  m/sec. 


z 

(m) 

A 

T 

(K) 

S  M(xlO^)  R 

(%)  (kg/m3)  (mm/hr) 

Z(xlGT3) 

(nm^/nr) 

R  - 

25  mm/ hr 

1500 

278.0 

278.0 

100.0 

1.24 

25.0 

0.33 

1000 

282.9 

282.8 

76.7 

1.23 

24.1 

0.24 

500 

287.8 

287.5 

60.5 

1.18 

23.2 

0.20 

0 

292.7 

292.1 

48.8 

1.12 

22.0 

0.18 

R  - 

50  mm/hr 

1500 

278.0 

278.0 

100.0 

2.24 

50.0 

0.92 

1000 

282.9 

282.8 

77.1 

2.22 

47.0 

0.49 

500 

287.8 

287.4 

61.6 

2.16 

45.2 

0.40 

0 

292.7 

291.8 

50.7 

2.07 

43.4 

0.36 

R  = 

75  mm/hr 

1500 

278.0 

278.0 

100.0 

3.28 

75.0 

1.77 

1000 

282.9 

282.8 

77.4 

3.25 

71.7 

0.76 

500 

287.8 

287.2 

62.7 

3.17 

69.1 

0.61 

0 

292.7 

291.5 

52.4 

3.06 

66.7 

0.55 

R  = 

100  mm/hr 

1500 

278.0 

278.0 

100.0 

4.03 

100.0 

2.51 

1000 

282.9 

282.7 

77.7 

3.99 

90.1 

0.94 

500 

287.8 

287.2 

63.4 

3.90 

87.7 

0.76 

0 

292.7 

291.3 

53.6 

3.78 

85.0 

0.69 
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Table  19.  Vertical  distribution  of  thermodynamic  and  . 

hydrometeoric  quantities  for  Case  D  in  a  sub¬ 
cloud  downdraft  of  15  m/sec. 


z 

(m) 

TD 

(K) 

T 

(K) 

S 

(7.) 

M(xl03)  R 

(kg/m3)  (mm/hr) 

Z(xl0"5) 

(nm^/m3) 

R  - 

25  mm/hr 

1500 

278.0 

278.0 

100.0 

1.24 

25.0 

0.33 

1000 

282.9 

282.9 

76.5 

1.23 

24.2 

0.25 

500 

287.8 

287.6 

59.9 

1.19 

23.4 

0.21 

0 

292.7 

292.3 

47.8 

1.14 

22.3 

0.19 

R  - 

50  mm/hr 

1500 

278.0 

278.0 

100.0 

2.24 

50.0 

0.92 

1000 

282.9 

282.8 

76.7 

2.22 

47.4 

0.53 

500 

287.8 

287.5 

60.6 

2.17 

45.6 

0.43 

0 

292.7 

292.1 

49.1 

2.09 

43.8 

0.38 

R  - 

75  mm/hr 

1500 

278.0 

278.0 

100.0 

3.28 

75.0 

1.77 

1000 

282.9 

282.8 

77.0 

3.25 

72.3 

0.82 

500 

287.8 

287.4 

61.4 

3.18 

69.5 

0.65 

0 

292.7 

291.9 

50.3 

3.08 

67.0 

0.57 

R  - 

100  mm/hr 

1500 

278.0 

278.0 

100.0 

4.03 

100.0 

2.51 

1000 

282.9 

282.8 

77.1 

3.99 

91.6 

1.02 

500 

287.8 

287.3 

61.9 

3.91 

88.2 

0.80 

0 

292.7 

291.7 

51.1 

3.80 

85.3 

0.71 

3.  For  a  given  downdraft,  the  thermodynamic  and  hydrometeoric 
quantities  in  Case  D  are  not  very  different  from  those  in  Cases  B  and  C 
for  small  initial  precipitation  rates  (25  and  50  nm/hr) ;  for  larger 
initial  precipitation  rates  (75  and  100  nm/hrX  the  differences  are 
greater.  This  is  due  to  the  smaller  number  of  larger  drops  at  small  R, 
as  indicated  in  Fig.  8. 

The  drop-size  distributions  resulting  from  Case  D  for  various 
constant  downdraft  velocities  are  given  in  Figs.  19  and  20  for  initial 
rainfall  rates  of  25  and  100  mm/hr,  respectively.  Comparing  the  large 
radius  ends  of  these  two  figures,  it  is  easily  seen  how  drastically  the 
concentration  of  large  drops  is  reduced  by  collisional  breakup.  This 
is  especially  true  in  Fig.  20.  The  sharp  variations  that  occur  near 
0.05  cm  and  0.14  cm  result  from  the  numerical  formulation  given  in 
Chapter  4,  in  which  the  loss  term  is  computed  differently  in  the  three 
segments  of  the  drop-size  range.  As  indicated  earlier,  these  sharp 
breaks  probably  do  not  occur  in  nature.  They  were  employed  to  permit 
comparisons  could  be  made  with  the  results  of  List  and  Gillespie  (1976) 
and  Gillespie  and  List  (1978/79). 

Figs.  21  and  22  compare  the  drop-size  distributions  in  Cases  C 
and  D  for  precipitation  rates  of  25  and  100  mm/hr,  respectively.  Fig. 

21  illustrates  how  collisional  breakup  depletes  the  large  drops  and 
distributes  the  fragments  into  the  smaller  radius  categories.  The 
concentration  of  small  drops  in  Fig.  21  approaches  the  same  value  for 
both  Case  C  and  Case  D  because  evaporation  proceeds  quickly  for  these 
drops,  while  drops  in  the  radius  range  from  0.03  cm  to  0.12  cm  sweep  out 
the  drops  of  radii  less  than  0.03  cm  due  to  collision-coalescence. 


except  for  Case  D 
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This  effect  is  more  dramatic  in  Fig.  22  where  in  Case  D  there  are  fewer 
drops  both  with  radii  less  than  0.05  cm  and  with  radii  greater  than 
0.15  cm  than  in  Case  C.  These  results  indicate  that  rainshafts  with 
weak  downdrafts  may  result  in  ground  base  distributions  having  fewer 
drops  at  both  the  smaller  and  the  larger  ends  of  the  distribution  when 
the  precipitation  rate  is  higher. 

Figs.  23  and  24  depict  the  evolution  of  drop-size  distributions 
in  Case  D  for  initial  precipitation  rates  of  25  and  100  mm/hr, 
respectively.  Collisional  breakup  proceeds  quickly  as  noted  by  the 
large  decrease  in  the  concentration  of  large  drops  in  the  first  500  m 
below  cloud  base.  The  fragments  are  distributed  in  all  drop  categories 
below  a  drop  radius  of  0.15  cm.  The  effect  of  evaporation  is  not  as 
strong  for  the  500-m  distribution  as  for  the  1000-  and  1500-m  distri¬ 
butions.  Downward  from  the  cloud  base,  the  concentration  of  larger 
drops  can  be  seen  to  decrease  so  that  the  number  of  small 
fragments  produced  also  decreases  and  the  effect  of  evaporation  becomes 
more  apparent.  Fig.  24  illustrates  the  same  effect,  except  more 
dramatically.  Also,  Fig.  24  shows  the  effect  of  collision-coalescence 
of  the  mid-size  drops  on  the  smaller  drops. 

When  Figs.  23  and  24  were  analyzed,  a  significant  difference  was 
noted  between  the  Case  D  results  and  the  results  reported  by  List 
and  Gillespie  (1976)  and  Gillespie  and  List  (1978/79) .  Replotted  in 
Fig.  25  are  the  results  of  Gillespie  and  List  (1978/79)  on  the  evolution 
of  the  drop-size  spectrum  through  collision-coalescence  and  collisional 
breakup  for  an  initial  precipitation  rate  of  25  mm/hr.  The  shape  of 
the  distribution  at  the  ground  (Z  ■  0  m)  in  Fig.  25  agrees  fairly 


Log  Drop  Concentration  (m  cm 


Radius  (cm) 


Fig.  23.  Evolution  of  drop-size  distribution  in  Case  D  for  an 
initial  precipitation  rate  of  25  mm/hr  and  a  constant 
downdraft  of  5  m/sec. 
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-  x*2km 


Fig.  25.  Evolution  with  altitude  of  drop-size 

distribution  through  collision-coalescence 
and  collisional  breakup  from  an  initial 
Marshall-Palmer  distribution  for  rainfall 
rate  of  25  mm/hr  (after  Gillespie  and  List, 
1978/79). 


closely  with  the  shape  of  the  1500-m  distribution  in  Fig.  23,  but  the 
drop  concentration  at  the  peak  of  the  distribution  is  higher  in  Fig.  25 
than  in  Fig.  23.  In  Fig.  26,  in  which  the  results  of  Gillespie  and 
List  are  plotted  against  those  of  Case  D  at  1000  m  and  1500  m  below 
cloud  base,  the  peaks  in  the  Gillespie  and  List  curves  cross 
the  Marshall-Palmer  distribution  at  a  drop  radius  of  about  0.075  cm, 
while  in  Case  D  the  crossover  occurs  at  about  r  *  0.14  cm.  Gillespie 
and  List  did  not  report  results  for  heavier  precipitation  rates,  al¬ 
though  in  their  earlier  work,  List  and  Gillespie  (1976)  provided  the 
distribution  resulting  from  an  initial  precipitation  rate  of  100  mm/hr. 
The  comparison  of  our  results  with  this  case  is  discussed  below. 

Fig.  27  is  a  time  evolution  through  coalescence  and  collisional 
breakup  of  an  initial  Marshall-Palmer  spectrum  for  a  rainfall  rate  of 


Concentration 


Fig.  26.  Drop-size  distributions  in  Case  D  compared  with 

analogous  computation  of  Gillespie  and  List  (1978/ 
79). 
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Fig.  27.  Time  evolution  through  coalescence  and 
breakup  of  a  raindrop  spectrum  from  an 
initial  Marshall-Palmer  distribution 
for  the  rainfall  rate  of  100  mm/hr 
(after  List  and  Gillespie,  1976) . 


100  mm/hr  as  reported  by  List  and  Gillespie  (1976) .  The  drop- 
concentration  characteristics  that  appear  in  this  figure  are  similar 
to  those  in  Fig.  25,  and  if  compared  against  Fig.  24,  would  reveal 
differences  from  Case  D  as  seen  above  for  a  25-tnm/hr  precipitation 
r»te.  k  possible  explanation  for  the  difference  is  in  the  formulation 
of  fragment  probability  resulting  from  collisional  breakup.  It  appears 
there  is  an  error  in  Eq.  (17)  of  List  and  Gillespie.  The  Gaussian 


distribution  fitted  to  the  observed  points  for  large  fragments  did  not 
have  a  negative  sign  in  the  exponent  -  an  error  that  this  research 
corrected  in  (40).  At  the  time  the  error  was  noted,  it  appeared  to  be 
restricted  to  the  journal  only  and  not  the  actual  List  and  Gillespie 
computer  program.  In  order  to  determine  the  effect  of  this  error 
on  the  final  drop-size  distribution  if,  in  fact.  List  and  Gillespie  had 
programmed  their  (17)  as  reported  in  the  literature,  a  test  case  was 
run  for  an  initial  precipitation  rate  of  100  mm/hr  and  a  constant 
downdraft  of  5  m/sec.  The  results  of  this  test  are  shown  in  Fig.  28. 
Since  List  and  Gillespie  (1976)  did  not  employ  rainshaft  downdraft 
velocity  in  their  research,  the  case  of  w^  *  5  m/sec  was  chosen  for 
comparison,  this  being  the  case  with  the  weakest  downdraft  in  our  com¬ 
putations.  The  results  in  Fig.  28  show  that  the  peak  concentration 
and  the  crossover  points  of  the  Marshall-Palmer  distribution  agree  quite 
well  between  the  test  case  (the  modified  Case  D)  and  List  and  Gillespie 
(1976).  The  only  major  disagreement  in  these  two  figures  is  the  drop 
radius  of  the  peak  concentration.  This  peak  occurs  at  r  -  0.03  cm  for 
List  and  Gillespie  and  r  *  0.05  cm  for  the  modified  Case  D.  This  is 
probably  due  to  the  numerical  formulation  used  by  List  and  Gillespie 
(1976). 

Since  there  was  strong  evidence  of  an  error  in  List  and 


Gillespie  (1976),  that  probably  was  carried  over  to  Gillespie  and 
List  (1978/79),  this  study  continued  to  the  final  step  of  adding 
aerodynamic  breakup  to  the  Case  D  microphysical  processes. 


Concentration 
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e.  Case  E:  Evaporation.  Collision-Coalescence,  Collisional  and 
Aerodynamic  Breakup 

Case  E  includes  the  Case  D  microphysical  processes  and  formula¬ 
tions  together  with  aerodynamic  breakup  formulations  of  Srivastava 
(1971)  and  is  the  last  case  considered  in  this  research. 

Vertical  distributions  of  thermodynamic  and  hydrometeoric 
quantities  resulting  from  Case  E  for  the  initial  precipitation  rates 
of  25,  50,  75,  and  100  mm/ hr  are  given  in  Tables  20,  21,  and  22, 
respectively,  for  constant  downdrafts  of  5,  10,  and  15  m/sec.  It  may 
be  recalled  that  in  the  comparison  of  Cases  B  and  C,  the  effects  of  aero¬ 
dynamic  breakup  were  found  to  be  insignificant  for  small  precipitation 
rates,  but  became  significant  with  larger  precipitation  rates  with 
their  higher  concentration  of  large  drops.  There  is  virtually  no 
difference  between  Case  D  and  the  Case  E  even  for  the  larger  precipi¬ 
tation  rates.  Collisional  breakup  depletes  the  concentration  of  large 
drops  to  the  point  where  aerodynamic  breakup  becomes  an  insignificant 
process . 

Drop-size  distributions  resulting  from  Case  E  are  given,  respec¬ 
tively,  in  Figs.  29,  30,  31,  and  32  for  initial  precipitation  rates  of 
25,  50,  75,  and  100  mm/hr.  A  comparison  of  Figs.  29  and  32  with  Figs. 

19  and  20  indicates  no  perceptible  difference  between  the  drop-size 
distributions  for  Cases  D  and  E,  thereby  showing,  again,  the  insignif¬ 
icance  of  aerodynamic  breakup  compared  to  collisional  breakup.  The 
distributions  in  Figs.  29-32  maintain  a  Marshall-Palmer  (exponential) 
slope  between  the  drop  radii  of  0.05  to  0.14  cm.  For  a  light  init¬ 
ial  precipitation  rate  (Fig.  28),  the  final  distributions  maintain  this 


Table  20.  Vertical  distribution  of  thermodynamic  and 

hydrometeorlc  quantities  in  a  subcloud  down- 
draft  of  5  m/sec. 


z  Tn  T  S  M(xl03)  R  Z(xl(T5) 

(m)  <K7  (K)  (%)  (kg/m3)  (tmn/hr)  (mm6/m3) 


R  * 

25  mm/hr 

1500 

278.0 

278.0 

100.0 

1.24 

25.0 

0.33 

1000 

282.9 

282.8 

77.4 

1.22 

24.0 

0.22 

500 

287.8 

287.3 

62.4 

1.18 

23.1 

0.19 

0 

292.7 

291.6 

51.8 

1.11 

22.0 

0.17 

R  » 

50  mm/hr 

1500 

278.0 

278.0 

100.0 

2.24 

50.0 

0.92 

1000 

282.9 

282.7 

78.1 

2.21 

46.7 

0.45 

500 

287.8 

287.0 

64.4 

2.15 

45.1 

0.38 

0 

292.7 

291.1 

55.1 

2.07 

43.5 

0.35 

R  ■ 

75  mm/hr 

1500 

278.0 

278.0 

100.0 

3.28 

75.0 

1.77 

1000 

282.9 

282.6 

78.8 

3.25 

71.2 

0.68 

500 

287.8 

286.8 

66.3 

3.17 

69.1 

0.58 

0 

292.7 

290.7 

58.2 

3.07 

67.2 

0.54 

R  - 

100  nm/hr 

1500 

278.0 

278.0 

100.0 

4.03 

100.0 

2.51 

1000 

282.9 

282.6 

79.2 

3.99 

90.1 

0.84 

500 

287.8 

286.6 

67.5 

3.91 

87.8 

0.72 

0 

292.7 

290.4 

60.2 

3.79 

85.7 

0.69 

Table  21.  Vertical  distributions  of  thermodynamic  and 
hydrometeoric  quantities  for  Case  E  in  a  sub¬ 
cloud  downdraft  of  10  m/sec. 


S  M(xl03) 
(7.)  (kg/m3) 

25  mm/hr 


R 

(mm/hr) 


Z(xl0“5) 

(nm^/m3) 


1500 

278.0 

278.0 

100.0 

1.24 

25.0 

0.33 

1000 

282.9 

282.8 

76.7 

1.23 

24.1 

0.24 

500 

287.8 

287.5 

60.5 

1.18 

23.2 

0.20 

0 

292.7 

292.1 

48.8 

1.12 

22.0 

0.18 

R  - 

50  mm/hr 

1500 

278.0 

278.0 

100.0 

2.24 

50.0 

0.92 

1000 

282.9 

282.8 

77.1 

2.22 

47.0 

0.50 

500 

287.8 

287.4 

61.6 

2.16 

45.2 

0.40 

0 

292.7 

291.8 

50.7 

2.07 

43.3 

0.36 

R  - 

75  mm/hr 

1500 

278.0 

278.0 

100.0 

3.28 

75.0 

1.77 

1000 

282.9 

282.8 

77.4 

3.25 

71.7 

0.76 

500 

287.8 

287.2 

62.7 

3.17 

69.0 

0.61 

0 

292.7 

291.5 

52.4 

3.06 

66.6 

0.55 

R  - 

100  msr/ihr 

1500 

278.0 

278.0 

100.0 

4.03 

100.0 

2.51 

1000 

282.9 

282.7 

77.7 

3.99 

90.8 

0.94 

500 

287.8 

287.1 

63.4 

3.90 

87.6 

0.75 

0 

292.7 

291.3 

53.6 

3.78 

84.8 

0.69 

Table  22.  Vertical  distributions  of  thermodynamic  and 

hydrometeorlc  quantities  for  Case  E  in  a  sub¬ 
cloud  downdraft  of  15  m/sec. 


z 

(m) 

td 

00 

T 

00 

S  M(xlO^)  R  Z(xl0“5) 

(7.)  (kg/m^)  (mm/hr)  (mm^/nP) 

R  - 

25  mm/ hr 

1500 

278.0 

278.0 

100.0 

1.24 

25.0 

0.33 

1000 

282.9 

282.9 

76.5 

1.23 

24.2 

0.25 

500 

287.8 

287.6 

59.9 

1.19 

23.4 

0.21 

0 

292.7 

292.3 

47.8 

1.14 

22.3 

0.18 

R  - 

50  mm/hr 

1500 

278.0 

278.0 

100.0 

2.24 

50.0 

0.92 

1000 

282.9 

282.8 

76.7 

2.22 

47.4 

0.53 

500 

287.8 

287.5 

60.6 

2.17 

45.5 

0.43 

0 

292.7 

292.1 

49.1 

2.09 

43.6 

0.38 

R  - 

75  mm/hr 

1500 

278.0 

278.0 

100.0 

3.28 

75.0 

1.77 

1000 

282.9 

282.8 

77.0 

3.25 

72.3 

0.82 

500 

287.8 

287.4 

61.4 

3.18 

69.4 

0.64 

0 

292.7 

291.9 

50.3 

3.08 

66.9 

0.57 

R  - 

100  mm/hr 

1500 

278.0 

278.0 

100.0 

4.03 

100.0 

2.51 

1000 

282.9 

282.8 

77.1 

3.99 

91.6 

1.02 

500 

287.8 

287.3 

61.9 

3.91 

88.0 

0.80 

0 

292.7 

291.7 

51.1 

3.80 

85.1 

0.71 

Log  Drop  Concentration  (m  cm 


Fig.  29.  Same  as  Fig.  7  except  for  Case  E 
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exponential  characteristic  down  to  a  radius  of  about  0.02  cm.  Young 
(1975)  and  Srivastava  (1978)  achieved  similar  results  when  they  stud¬ 
ied  the  time-dependent  interaction  of  condensation,  coalescence,  and 
collisional  and  aerodynamic  breakup  in  the  initiation  of  warm  rain. 
Both  researchers  used  the  aerodynamic  breakup  formulated  by  Srivastava 
(1971)  and  the  collision- induced  breakup  of  Brazier-Smith  et  aJL.  (1972) 
in  their  studies.  Young  found  that  not  only  did  collision  breakup 
dominate  over  aerodynamic  breakup,  but  the  resulting  steady-state  dis¬ 
tribution  was  exponential  and  provided  fair  agreement  with  average 
observed  drop  spectra.  Srivastava  (1978)  appeared  to  find  similar 
results  in  his  theoretical  study. 
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CHAPTER  VI 

COMPARISON  WITH  OBSERVED  RAINDROP  SPECTRA 

We  have  computed  steady-state  raindrop-size  distributions  for 
several  constant  downdraft  velocities  and  cloud-base  precipitation 
rates,  but  how  well  do  these  calculated  distributions  agree  with 
observed  spectra?  In  this  chapter  we  will  compare  our  computed  re¬ 
sults  with  measured  spectra  from  both  warm  rain  and  also  thunder¬ 
storm  (cold  rain)  situations. 

a.  Warm  Raindrop  Measurements 

Austin  and  Geotis  (1979)  reported  several  sets  of  drop-size 
measurements  made  on  board  a  ship  during  the  Global  Atmospheric 
Research  Program's  (GARP)  Atlantic  Tropical  Experiment  (GATE).  They 
used  a  Joss  disdrometer  which  was  manually  activated  during  moderate 
and  heavy  showers.  The  Joss  disdrometer  has  the  capability  to  record 
drop  sizes  in  the  radius  range  0.02  to  about  0.25  cm,  but  due  to  noisy 
conditions  aboard  ship,  most  drops  <  0.10  cm  in  radius  were  suppressed 
in  the  measurements.  For  this  reason  very  large  correction  factors 
had  to  be  added  to  the  small  drop  measurements.  For  a  full  explana¬ 
tion  of  the  problems  and  solutions  associated  with  these  measurements, 
see  Austin  and  Geotis  (1979). 

Fig.  33  is  plot  of  an  Austin-Geotis  drop  distribution  for 
R  *  34.4  mm/hr  together  with  Case  E  distributions  for  initial  precipi¬ 
tation  rates  of  25  and  50  mm/hr.  The  observed  distribution  has  more 


drops  of  the  small  radii  and  fewer  at  the  large  radii  than  either  of 


Radius  (cm) 

Fig.  33.  Comparison  of  drop-size  distributions  computed  in  Case  E 
with  those  observed  by  Austin  and  Geotis  (1979)  in 
convective  rain.  Computed  distributions  refer  to  1500  m 
below  cloud  base,  and  are  for  initial  precipitation  rates 
of  25  and  50  mm/hr  and  a  downdraft  velocity  of  5  m/sec. 
Observed  distributions  are  for  a  precipitation  rate  of 
34.4  mm/hr  at  the  surface. 


the  computed  distributions.  In  Fig.  34,  a  comparison  of  distributions 
resulting  from  heavier  rain  shows  that  the  observed  distribution  has 
the  same  general  shape,  but  overall  a  higher  concentration  of  drops 
than  the  computed  spectra.  Both  Figs.  33  and  34  are  for  computed 
distributions  resulting  from  a  5  m/sec  constant  downdraft.  Fig.  35 
shows  that  the  Case  E  spectrum  for  wQ  =  15  m/sec  is  in  closer  agree¬ 
ment  with  Austin  and  Geotis  than  the  two  distributions  resulting  from 
weaker  downdrafts.  This  is  especially  true  in  the  radius  range  from 
0.14  to  0.18  cm.  This  result  should  be  expected  since  observations 
were  taken  in  moderate  to  heavy  rainshafts  where  the  downdraft  velo¬ 
city  could  be  stronger. 

In  order  to  insure  that  their  instrument  corrections  were  valid, 
Austin  and  Geotis  compared  their  observed  spectra  with  distributions 
observed  in  other  maritime  regions.  The  spectra  from  Majuro  Atoll  in 
the  Marshall  Islands,  observed  by  Mueller  and  Sims  (1967\  and  the 
Panama  drop-size  distributions  measurements  made  by  Geotis  i,3^68)  wstf'**-: 
selected  for  comparison.  Mueller  and  Sims  used  a  raindrop  camera, 
while  Geotis  (1968)  used  a  Joss  disdrometer.  Fig.  36  is  a  plot  of 
averaged  drop-size  distributions  measured  by  Austin  and  Geotis  (1979) , 
Mueller  and  Sims  (1967),  and  Geotis  (1968)  together  with  a  computed 
distribution  at  1500  m  below  cloud  base  for  Case  E  and  an  initial 
precipitation  rate  of  25  mm/hr.  The  computed  distribution  has  a 
higher  drop  concentration  across  the  entire  radius  spectrum  than  any 
of  the  averaged  observed  distributions.  This  trend  is  less  noticeable 
with  higher  precipitation  rates  as  shown  in  Fig.  37.  Here,  the  cal¬ 
culated  drop  concentration  agrees  much  better  with  average  observed 
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Fig.  36.  Averaged  observed  drop-size  distributions  from  Mueller 
and  Sims  (1967),  Geotis  (1968),  and  Austin  and  Geotis 
(1979)  together  with  a  computed  distribution  for  Case 
E.  The  initial  precipitation  rate  is  25  mm/hr. 


110 


spectra.  Note  In  both  Figs.  36  and  37  that  the  Geotls  (1968)  distri¬ 
butions  have  a  much  higher  drop  concentration  for  small  radius  drops. 
This  is  in  the  minimum  radius  range  of  drop  resolution  for  the  Joss 
disdrometer  and,  based  on  acoustical  problems  in  this  radius  range, 
should  probably  not  be  highly  regarded.  Two  things  are  clear,  how¬ 
ever,  from  Figs.  33  through  37.  First,  the  computed  distributions 
for  R  *  25  mn/hr  have  drop  concentrations  considerably  higher  than 
observed  spectra  for  the  same  precipitation  rate,  while  Figs.  34,  35, 
and  37  indicate  that  for  higher  precipitation  rates  this  trend  be¬ 
comes  less  evident.  Second,  all  distributions  are  very  steep,  with 
numerous  drops  in  the  range  0.05  cm  <_  r  <_  0.10  cm  and  relatively  few 
drops  with  r  >  0.15  cm.  This  last  conclusion  is  important  since  it 
shows  that,  in  warm  rain,  collislonal  drop  breakup  is  the  driving 
mechanism  in  shaping  the  final  drop-size  distribution. 

b.  Thunderstorm  Raindrop  Measurements 

Thunderstorm  raindrop-size  distributions  were  measured  by  Dingle 
and  Hardy  (1962)  using  a  photoelectric  raindrop  spectrometer.  They 
took  several  measurements  from  storms  in  the  fall  of  1959  in  the  Ann 
Arbor,  Michigan,  area  that  would  strongly  indicate  frozen  hydrometeors 
in  the  rainshaft.  Fig.  38  depicts  an  averaged  Dingle  and  Hardy  dis¬ 
tribution  (average  of  seven  1-mln  distributions)  for  R  ■  24.2  mm/hr 
plotted  against  the  distributions  in  Case  G  for  R  ■  25  nm/hr  and 
Wp  *  5  and  15  m/sec.  As  in  the  warm  rain  comparison  (Fig.  33),  the 
observed  distribution  has  a  lower  concentration  of  large  drops  when 
compared  to  the  calculated  distribution,  although  this  difference  is 
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Fig.  38.  Drop-size  distributions  computed  in  Case  E  for  an 

initial  precipitation  rate  of  25  mm/hr  compared  with 
an  averaged  drop-size  distribution  observed  by  Dingle 
and  Hardy  (1962)  in  thunderstorm  rain  of  24.2  mm/hr 
intensity. 


less  in  Fig.  38  than  in  Fig.  33.  The  difference  between  observed  and 
calculated  drop  concentrations  increases  in  Fig.  38  toward  smaller 
drop-sizes.  Note  the  increase  in  drop  concentration  (Fig.  38)  which 
occurs  at  about  r* 0.17  cm.  As  Indicated  earlier,  the  Dingle  and 
Hardy  distributions  are  averages  which  means  this  has  to  be  a  persis¬ 
tent  phenomenon.  A  similar  trend  is  apparent  in  Fig.  33  at  about 
r  ■  0.15  cm.  The  comparison  of  Figs.  33  and  38  illustrates  a  greater 
effect  o£  drop  breakup  on  warm  rain  than  on  thunderstorm  rain. 

Fig.  39  shows  a  comparison  of  averaged  Dingle  and  Hardy  distributions 
(average  of  three  1-mln  distributions)  for  R  ■  65.5  mm/hr  plotted 
against  Case  E  results  for  R  ■  75  mm/ hr.  The  trend  in  Fig.  39  is 
the  same  as  in  Fig.  38  in  that  the  observed  drop  concentrations  are 
lower  than  Case  E  values  until  about  r  -  0.18  cm,  where  the 
slope  of  the  observed  distribution  changes  drastically.  Having  no 
detail  is  available  on  the  drop-size  resolution  of  the  equipment  or 
on  any  problems  associated  with  the  measurements  of  Dingle  and  Hardy 
it  will  not  be  proper  to  speculate  on  the  possible  cause  of  this  drop 
concentration  increase. 

In  summary,  in  the  radius  range  from  0.03  to  about  0.19  cm,  the 
calculated  drop-size  distributions  appear  to  agree  with  the  observed 
average  distributions  in  shape,  but  disagree  in  drop  concentration  in 
the  small  drop-size  portions  of  the  spectra. 
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CHAPTER  VII 


SUMMARY  AND  CONCLUSIONS 

By  assuming  steady  state  and  adiabatic  (nonentraining)  conditions 
in  a  constant  downdraft  below  a  cloud  with  raindrops  at  the  cloud  base 
following  a  Marshall-Palmer  distribution,  we  have  computed  the  tem¬ 
perature,  relative  humidity,  and  drop-size  distribution  in  the  down- 
draft  as  it  progresses  downward.  Modification  in  the  drop-size 
distribution  has  been  considered  as  resulting  from  the  microphysical 
processes  of  evaporation,  collision-coalescence,  and  collisional  and 
aerodynamic  breakup.  The  drops  resolved  in  the  computation  cover  the 
range  from  0.004  cm  to  0.4  cm  in  radius;  the  drops  with  radii  less  than 
0.004  cm  are  collectively  treated  as  "cloud  droplets."  Drops  with 
radii  greater  than  0.4  cm  are  assumed  to  break  up  spontaneously  and 
instantly.  We  also  have  compared  our  calculated  drop-size  distribu¬ 
tions  with  several  measured  on  the  ground  and  on  ships . 

Evaporation  by  itself  (Case  A)  affects  primarily  drops  with  radii 
less  than  0.015  cm.  The  evaporation  time  constants  for  these  smaller 
drops  are  so  short  that  they  do  not  survive  a  height  interval  Az  ■  25  m 
once  the  relative  humidity  decreases  to  a  value  of  about  60%.  The 
concentration  of  the  small  drops  is  replenished  by  the  evaporation 
of  larger  drops,  but  since  the  evaporation  rate  is  less  for  larger 
drops,  fewer  drops  enter  and  leave  the  small  radius  categories.  Thus 
our  final  drop-size  distribution  peak  at  radii  from  0.015  to  0.020 
cm  depending  on  the  strength  of  the  constant  downdraft.  The 


peak  drop  concentration  shifts  to  smaller  radius  categories  for  strong-- 
er  downdrafts  since  the  faster  moving  drops  have  less  time  to  evapor¬ 
ate.  The  rainshaft  relative  humidity  is  minimum  and  temperature 
maximum  for  small  precipitation  rates  and  strong  downdrafts.  The 
minimum  relative  humidity  computed  was  47.8%  at  1500  m  below  cloud 
base  for  R  ■  25  mm/hr  and  wQ  ■  15  m/sec.  Evaporation  does  not 
affect  rainfall  rate  and  radar  reflectivity  to  a  great  extent,  the 
magnitude  of  the  decrease  was  less  than  lu%  in  the  most  severe  case. 

In  Case  B,  when  collision-coalescence  is  considered  together  with 
evaporation,  several  interesting  changes  occur.  The  precipitation 
rate  increases  dramatically  as  the  large  drops  with  their  higher 
terminal  velocities  sweep  out  the  small  drops,  while  the  radar  reflec¬ 
tivity,  which  depends  on  the  concentration  of  larger  drops  increases 
significantly.  The  maximum  Increases  occur  with  high  Initial  rain¬ 
fall  rates,  which,  according  to  the  Marshall-Palmer  distribution,  pro¬ 
vide  higher  initial  concentrations  of  large  drops,  and  with  weak 
downdrafts  permitting  the  larger  drops  a  longer  time  to  sweep  out 
the  smaller  drops.  Rainshaft  relative  humidity  is  lower  and  the 
temperature  higher  for  Case  B  than  Case  A  because  the  larger  drops 
which  are  poor  evaporators  deplete  the  more  efficiently  evaporating 
small  drops.  For  high  initial  rainfall  rates,  the  concentration  of 
drops  with  radii  less  than  0.075  cm  is  severely  depleted  by  collision- 
coalescence. 

In  Case  C,  when  aerodynamic  breakup  is  added  to  evaporation  and 
collision-coalescence,  one  notes  little  difference  from  Case  B  for 
small  initial  rainfall  rates.  The  effect  of  aerodynamic  breakup 


becomes  more  apparent  for  higher  initial  precipitation  rates  because 
of  the  higher  concentration  of  large  drops.  Subcloud  relative  humidity 
is  lower  and  temperature  higher  for  Case  C  when  compared  to  Case  B, 
although  this  effect  is  minor  for  small  initial  rainfall  rates. 

The  probability  of  aerodynamic  breakup  becomes  very  small  (<0. 01/sec) 
for  drops  of  radii  less  than  0.3  cm.  As  shown  earlier  in  Chapter 
VI,  the  fragment  population  resulting  from  aerodynamic  breakup  has  a 
maximum  in  our  spectrum  at  a  drop  radius  of  0.055  cm.  For  large 
initial  precipitation  rates,  the  drop-size  distribution  in  Case  C 
has  a  much  larger  concentration  of  small  drops  than  in  Case  B  because 
of  the  breakup  of  the  largest  drops. 

Case  D  combines  collisional  breakup  with  the  processes  of  eva¬ 
poration  and  collision-coalescence;  aerodynamic  breakup  is  not  in¬ 
cluded  in  this  case.  Collisional  breakup  is  a  much  more  efficient 
process  for  depleting  large  drops  than  is  aerodynamic  breakup.  As  a 
result,  the  subcloud  relative  humidity  is  higher  and  temperature 
lower  in  Case  D  than  in  Case  C  for  weak  downdrafts.  For  strong 
downdrafts,  this  effect  is  not  as  great  because  the  higher  concen¬ 
tration  of  small  drops  is  given  less  time  to  evaporate.  Radar 
reflectivity  as  well  as  rainfall  rate  in  the  subcloud  layers  is 
drastically  reduced  in  Case  D  due  to  rapid  breakup  of  the  large 
radius  drops.  Drop-size  distributions  in  this  case  exhibit  a  rapid 
decrease  in  drop  concentration  beginning  with  drops  of  radius 
0.15  cm  such  that  very  few  drops  with  radii  greater  than  0.20  cm 


exist  in  the  final  distribution. 


117 


In  Case  E,  when  aerodynamic  breakup  is  added  Co  Che  Case  D 
processes,  virCually  no  change  is  noCed.  This  is  due  Co  Che  rapid 
depleCion  of  drops  wich  large  radii  by  collisional  breakup.  Aero¬ 
dynamic  breakup  is  probably  a  more  imporcanc  microphysical  process 
in  cold  rain  resulcing  from  sCrong  convecCive  sCorms  and  jusC  be- 
neach  Che  melcing  level  where  large  drops  resulcing  from  Che  melting 
of  ice  pelleCs  and  hail  can  become  unsCable  and  break  up. 

The  compuCed  drop-size  discribucions  obcained  in  Case  E  compare 
favorably  wich  maricime  specCra  measured  by  AusCin  and  GeoCis  (1979), 
Mueller  and  Sims  (1967),  and  Geotis  (1968).  The  general  shape  of  Che 
compuCed  distributions,  although  quite  similar  to  measured  warm  rain 
spectra,  generally  shows  a  higher  concentration  of  drops  with  radii 
greater  than  0.10  cm  than  observed  with  small  precipitation.  How¬ 
ever,  for  higher  precipitation  rates  the  computed  distributions 
provide  an  excellent  fit  to  observed  spectra  in  the  radius  range 
0.05  -  0.20  cm,  which  was  the  range  of  the  measuring  disdrometers. 
Case  E  drop-size  distributions  also  were  compared  with  the  measured 
cold  rain  distributions  of  Dingle  and  Hardy  (1962) .  Computed  dis¬ 
tributions  displayed  higher  concentrations  of  small  drops  in  the 
radius  range  of  0.05  -  0.15  cm  than  in  measured  cold  rain  spectra. 
Although  the  measurement  of  drop  concentration  was  limited  at  the 
large  radius  end  (r>0.20  cm),  indications  are  that  the  observed  con¬ 
centration  of  drops  with  radii  greater  than  0.20  cm  would  be  higher 
than  the  computed  distributions. 
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In  conclusion,  the  study  reported  herein  throws  considerable 
light  on  the  role  of  the  microphysical  processes  in  determining  the 
thermodynamic  and  microphysical  characteristics  of  the  warm  rain 
accompanied  by  a  convective  downdraft.  The  results  are  encouraging 
because  of  the  qualitative  agreement  of  the  drop-size  distributions 
observed  in  nature  with  those  computed.  However,  there  should  be 
no  doubt  that  this  agreement  is  somewhat  fortuitous,  in  view  of  the 
extreme  simplicity  of  the  model  vis-a-vis  the  possibility  that  the 
details  of  the  basic  model  (such  as  the  Marshall-Palmer  drop-size 
distribution  at  the  cloud  base,  assumed  at  a  height  of  1500  m  above 
the  ground)  are  different  from  the  conditions  under  which  the 
measured  drop-size  distributions  had  been  obtained.  Obviously  fur¬ 
ther  research  is  required,  with  an  improved  model  incorporating 
realistic  initial  and  environmental  conditions. 


CHAPTER  VIII 


RECOMMENDATIONS 

In  developing  this  study,  several  simplifying  assumptions  were 
made,  many  of  which  were  invoked  so  we  could  compae  our  results  with 
those  of  previous  studies  incorporating  the  same  assumptions.  In 
addition,  our  goal  was  to  explore  the  role  of  the  fundamental  pro¬ 
cesses  rather  than  reproduce  results  in  natural  situations.  For 
closer  comparison  with  observations,  the  present  computations  need 
several  improvements  as  indicated  below. 

1.  Smoother  coalescence  efficiencies  utilizing  the  latest 
experimental  results  should  help  eliminate  the  sharp  breaks  which 
occur  at  radii  of  0.05  cm  and  0.15  cm.  Also  the  collection  kernel 
should  incorporate  more  realistic  collision  efficiencies  than  geome¬ 
tric  sweepout. 

2.  In  this  study,  it  was  assumed  that  breakup  was  constant  below 
the  radius  of  0.025  cm,  since  the  McTaggart-Cowan  and  List  (1975)  experi¬ 
ments  had  a  lower-size  limitation  of  r  *  0.025  cm.  The  results  of 
Brazier-Smith  et  al.  (1972)  included  drop  interactions  down  to  drop 
radii  of  about  0.015  cm.  Incorporating  the  experimental  data  of 
Brazier-Smith  et  al.  into  those  of  McTaggart-Cowan  and  List  (1975) 
could  make  the  resulting  drop-size  distributions  more  representative 

at  the  smaller  drop  radii. 

3.  Variable  subcloud  density  and  raindrop  terminal  velocity  in 
future  computations  should  improve  the  realism  of  the  computed  drop- 
size  spectra.  The  error  in  assuming  these  two  quantities  constant  is 

S 
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not  too  large  for  a  shallow  subcloud  region,  however. 

4.  The  assumptions  of  a  constant  subcloud  downdraft  and  of 
steady-state  are  unrealistic.  The  former  can  be  changed  easily 
for  better  realism.  The  latter  would  require  extensive  changes  in 
computational  procedure.  Introduction  of  the  effects  of  turbulence 
would  be  almost  prohibitively  difficult. 

5.  Cloud  base  drop-size  distributions  other  than  Marshall-Palmer 
should  be  investigated.  Specifically,  cloud  base  spectra  measured 
during  GATE  and  reported  in  the  literature  should  be  used  as  the 
initial  distribution.  This  will  enable  one  to  see  how  well  the 
resulting  steady  state  distributions  agree  with  those  measured  at  the 
surface  during  that  experiment . 
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APPENDIX  A 

LIST  OF  SYMBOLS 

Empirical  constant  (■  62.3) 

-3 

Liquid  water  content  in  precipitation  (kg  m  ) 
Empirical  constant  (*7) 

Specific  heat  of  air  at  constant  pressure 
Raindrop  diameter  (cm)  in  the  i-th  category 

2  -1 

Diffusivity  of  water  vapor  in  air  (m  sec  ) 

Power  law  constant  in  equation  (39) 

Logarithmic  difference  in  radii  between  any  two 
adjacent  size  classes 

Vapor  pressure  (Pa) 

Saturation  vapor  pressure  (Pa) 

Collection  efficiency  for  a  drop  of  radius  i  colliding 
with  a  drop  of  radius  1 

Collision  efficiency  for  drops  of  radius  i  and  1 
Coalescence  efficiency  for  drops  of  radius  i  and  1 
Total  number  of  small  fragments 

Dimensionless  factor  which  is  a  function  of  Reynolds 
Number 

2 

Acceleration  due  to  gravity  (9.8  m/sec  ) 

Size  category  of  drops 

Integer  subscript  Indicating  height  (jAz) 

Unit  vector  in  the  z  direction 
Size  category  of  drops 
Thermal  conductivity  of  air 


Adjustable  scale  factor  that  fixes  the  data  points  at 
integer  values  of  k 

The  class  index  of  drops  which  coalesce  with  drops 
between  class  kh  and  class  k  to  form  drops  of  class 
k+1 

The  class  index  of  drops  whose  mass  is  half  that  of 
class  k 

Latent  heat  of  vaporization 

Mass  of  a  single  drop  in  radius  class  i 

Number  of  drops  per  cubic  meter  in  radius  class  i 
(m*3) 


Number  of  drops  of  radius  r  per  cubic  meter  per  unit 
radius  class  interval  (m-3cm-3-) 

Total  number  of  drops  per  kilogram  of  air  (kg 

Number  of  drops  per  kilogram  of  air  (kg 

Number  of  drops  of  radius  r  per  kilogram  of  air  per 
unit  radius  interval  r 

Atmospheric  pressure  (Pa) 

Probability  of  a  drop  of  radius  r  to  be  formed  by 
the  collisional  breakup  of  a  large  (L)  and  a  small 
(s)  drop  per  unit  time  (sec-3-) 

Probability  of  a  drop  of  radius  r .  breaking  up  due  to 
instability  per  unit  time  (sec-3-)1 

Source  and  sink  term  in  the  continuity  equation  for 
drop  concentration 

Probability  by  which  a  breaking  drop  of  radius  r^  yields 
a  drop  of  radius  r^  per  unit  time  (sec-3-) 

Drop  radius  (cm) 

The  drop  radius  in  the  i-th  category  (cm) 

Precipitation  rate  (nrn/hr) 

Specific  gas  constant  for  dry  air 
Reynolds  number 
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W±(L,s) 


Specific  gas  constant  for  water  vapor 

Equivalent  thickness  of  the  transition  shell  outside 
the  raindrop 

Relative  humidity  (percent) 

Time  (sec) 

Environmental  temperature  (K) 

Mean  temperature  in  a  25-m  layer  (K) 

Dry  adiabatic  Temperature  (K) 

Velocity  of  air  (m/sec) 

Velocity  of  a  liquid  drop  of  radius  r  relative  to 
the  air  (m/sec) 

Terminal  speed  of  drops  (m/sec) 

Steady-state  downdraft  speed  (m/sec) 

Rate  of  formation  of  drops  of  radius  r^  by  colli- 

sional  breakup  of  drops  with  radii  r  and  r 

L  S 

Height  above  ground  (m) 

6  3 

Radar  reflectivity  (mm  /m  ) 

Radius  interval  of  the  i-th  category  (cm) 

Difference  between  drop  class  k  and  drop  class  kc 

rk+l/rk 

_3 

Dry  adiabatic  lapse  rate  (9.81  x  10  k/m) 

3 

Density  of  liquid  water  (kg/m  ) 

Thickness  of  a  layer  (25  m) 

Potential  temperature  (K) 

Exponent  in  Marshall-Palmer  Equation  (40) 

Kinematic  viscosity 

3 

Saturated  vapor  density  at  raindrop  surface  (kg/m  ) 


3 

Environmental  vapor  density  (kg/m  ) 

3 

Density  of  dry  air  (kg/m  ) 

3 

Liquid  water  content  in  air  (kg/m  ) 

Water  vapor  mixing  ratio 

1  dpd 

Defined  as - -3 — 

p  .  dz 

^  -1  -1 
Specific  entropy  (jkg  K  ) 

Collection  kernel  for  drops  of  radius  r^  colliding 
with  drops  of  radius  r^  (fflVsec) 
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APPENBIX  B 

ANALYSIS  OF  THE  ERROR  DUE  TO  THE  ASSUMPTION  OF  CONSTANT 
DENSITY  AND  CONSTANT  TERMINAL  SPEED 


Density 

The  advection  term  of  (12)  is 


dN 
_ r 

dz  EV 


-7-  (p  N  Vm) 
dz  r  T 


(B.l) 


Assuming  to  be  constant  in  the  subcloud  layer  and  setting 
a(2)  a  -  —  one  can  show  that 


(InN  )  -  - 
dz  r  »7 

w_+V_  . 
D  T* 


(B.2) 


which  is  solved  numerically  as 


"k.j-fi  ■  Vj 


w  +V 
D  T,± 


o(z)Az 


(B.  3) 


If  p  is  assumed  constant  over  Az,  a(z)  +  0,  and  N  .  .  =  M  . .  In 

k. ,  j  +1  k ,  j 

order  to  estimate  the  error  in  this  assumption,  (B.3)  is  solved 
for  large  raindrops  (r  *  0.33  cm)  by  using 

■  9.2  m/sec,  Az  *  1500  m,  a(z)  *  10  \  and  w^  *  5  m/sec. 

The  result  is  N^  ■  0.903  .  Assuming  p  to  be  a  constant 

over  Az  ■  1500  m  results  in  a  9.7%  error  in  N,  ...  . 

k,j+l 
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Terminal  Velocity 

Using  (B.l)  and  assuming  p  constant  in  the  subcloud  layer  we 

have 


m  .  «  Az  _ T 

k,j+l  k,j  SXP  wd+Vt  dz  .  (B.4) 

The  summer  atmosphere  results  of  Beard  (1976)  give  the  change  in 

dV  -1 

terminal  velocity  with  height,  _ T  ,  as  approximately  0.0005  sec 

dz 

for  drops  of  0.3  cm  radius.  With  Az  =  1500  m,  wQ  =  5  m/sec 

and  V„  =  9.2  m/sec, (B.4)  is  solved  as  N,  ...  =  0.947  N.  .. 

T  k,j+l  k,j 

The  error  incurred  by  assuming  a  constant  VT  in  the  subcloud  layers  is 

about  5.3%. 
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APPENDIX  C 

EXTENSION  OF  KINZER  AND  GUNN’S 
MASS  EVAPORATION  RESULTS 

The  rate  of  mass  evaporation  from  liquid  water  drops  in  motion 
relative  to  their  environment  has  been  explored  both  theoretically 
and  experimentally  by  Kinzer  and  Gunn  (1951)  .  Their  experimental 
results  are  expressed  as 

F 

-  =  [4irr(l  +  -f)}  [D(pa  -  ob)  ]  (C.l) 

and  the  two  factors  on  the  right-hand  side  are  shown  in  Table  Cl.  In 
(C.l),  s  is  the  equivalent  thickness  of  the  transition  shell  outside 
the  drop,  F  is  a  dimensionless  factor  which  is  a  function  of  Reynolds 
number  Re,  D  is  the  diffusivity  of  water  vapor  in  air,  Pa  the  saturated 
vapor  density  at  the  surface  of  the  drop,  and  is  the  vapor  density 
of  the  environment. 

This  research  deals  with  raindrops  of  radii  ranging  from  0.004  cm 
to  0.40  cm,  which  requires  an  extension  of  the  Kinzer  and  Gunn  data. 
Byers  (1965)  showed  that  the  mass  evaporation  rate  for  drops  in  motion 
could  be  written, 

”  *  [4irr(l  +  0.22  F  Re1/2)]  ,  (C.2) 

2  2  1/2 

where  a'  -  L  /KR  T  ,  b'  -  R  T/e  D,  and  (1  +  0.22  F  Re  '  )  is  the  venti- 
V  vs 

lation  factor.  Eq.  (C.2)  was  used  to  expand  the  evaporation  data 
since  (C.l)  required  a  knowledge  of  the  equivalent  thickness  of  the 
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transition  shell  outside  the  drop,  which  was  not  available. 


Table  C.  1.  Evaporation  of  freely  falling  water  drops  (after  Kinzer 
and  Gunn,  1951) . 


Table  C. 1A 


Table  C  .IB 


F 

4ur(l  +  — )* 
s 


D(p 

a 


oh)* 


Temperature  (C)  „  Temperature  (C) 

Radius  S  _ _ _ 


0 

10 

30 

40 

(7.) 

0 

10 

30 

40 

0.005 

0.086  0.082  0.079  0.076  0.073 

10 

0.61 

0.98 

1.47 

2.06 

2.68 

0.010 

0.29 

0.29 

0.29 

0.28 

0.28 

20 

0.54 

0.87 

1.29 

1.79 

2.36 

0.015 

0.49 

0.48 

0.48 

0.47 

0.47 

30 

0.48 

0.76 

1.12 

1.55 

2.05 

0.020 

0.73 

0.72 

0.71 

0.70 

0.69 

40 

0.41 

0.65 

0.95 

1.32 

1.75 

0.025 

1.01 

0.99 

0.97 

0.96 

0.94 

50 

0.34 

0.54 

0.78 

1.09 

1.45 

0.030 

1.31 

1.29 

1.27 

1.25 

1.24 

60 

0.27 

0.43 

0.63 

0.86 

1.15 

0.035 

1.66 

1.63 

1.61 

1.58 

1.55 

70 

0.20 

0.32 

0.46 

0.64 

0.85 

0.040 

2.03 

2.00 

1.97 

1.94 

1.91 

80 

0.14 

0.21 

0.31 

0.42 

0.56 

0.045 

2.5 

2.4 

2.4 

2.3 

2.3 

90 

0.07 

0.11 

0.16 

0.21 

0.28 

0.050 

2.9 

2.8 

2.8 

2.7 

2.7 

100 

0 

0 

0 

0 

0 

0.06 

3.9 

3.8 

3.7 

3.6 

3.6 

0.07 

4.9 

4.8 

4.7 

4.6 

4.5 

0.08 

6.0 

5.9 

5.8 

5.7 

5.6 

*See 

p.  133  for  explanation. 

0.09 

7.3 

7.2 

7.0 

6.9 

6.8 

0.10 

8.8 

8.5 

8.3 

8.1 

8.0 

0.11 

10.5 

10.1 

9.9 

9.6 

9.4 

0.12 

12.4 

12.0 

11.7 

11.3 

11.0 

0.13 

14.7 

14.2 

13.8 

13.3 

12.8 

0.14 

17.2 

16.6 

16.0 

15.4 

14.9 

0.15 

20.1 

19.3 

18.5 

17.8 

17.2 

0.16 

23 

22 

21 

21 

20 

0.17 

27 

26 

25 

24 

23 

0.18 

31 

30 

28 

27 

26 

o.:,9 

35 

34 

32 

31 

29 

0.20 

39 

38 

36 

35 

33 

0.21 

43 

42 

40 

39 

37 
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a.  Method  of  Computation 

The  first  step  was  to  create  Table  C.  2  for  the  same  radii  range 
as  in  Table  C  1.  This  was  done  by  first  calculating  Table  C.2B  in 
terms  of  the  second  bracketed  factor  in  (C.2)  with  temperature- 
dependent  values  of  latent  heat  of  vaporization,  thermal  conductivity, 
saturation  vapor  pressure,  and  diffusivity  taken  from  the  Smithsonian 
Meteorological  Tables  (List,  1949). 

Next,  Table  C. 2A,  for  the  radii  range  0.005  to  0.21  cm,  was  com¬ 
puted  by  using  the  following  process: 

1.  Ten  mass  evaporation  tables,  one  for  each  relative  humidity, 
were  constructed  by  multiplying  the  values  in  Table  C.  1A  by  those  in 
Table  C.1B  for  the  same  temperature;  values  in  the  table  represented 
(-  dM/dt)c^,  that  is,  the  mass  evaporation  rate  represented  by  Table 
C.l. 

2.  Ten  tables  similar  to  Table  C.2A,  one  for  each  relative 
humidity  value,  were  calculated  using 

4„r(l  +  0.22  F  Re1/2)  -  <C.3> 

3.  An  arithmetic  average  of  the  ten  tables  constructed  in  step 
(2)  is  shown  in  Table  C.  2A  for  drop  radii  0.005  to  0.21  cm. 

Table  C.2A  was  extended  to  radius  values  below  0.005  cm  by 
assuming  F  *  1,  which  is  consistent  with  the  values  given  by  Kinzer  and 
Gunn  (1951) .  The  Reynolds  number  data  required  for  these  calculations 
were  taken  from  List  (1949).  Table  C.  2A  was  extended  from  0.21-  to 
0.40-cm  radius  in  two  steps  as  follows. 
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Table  C.2.  Computed  evaporation  of  freely  falling  water  drops. 


Table  C.2A 


Table  C.2B 


4irr  (1+0,22  FRe1/2) 


S-l 

a'+b' 


Temperature  (C)  Temperature  (C) 


(cm) 

0 

10 

20 

30 

40 

(%) 

0 

10 

20 

30 

40 

0.004 

0.050 

0.038 

0.029 

0.023 

0.017 

10 

0.95 

1.91 

3.59 

6.28 

10.20 

0.005 

0.055 

0.042 

0.028 

0.024 

0.019 

20 

0.85 

1.70 

3.19 

5.58 

9.08 

0.010 

0.19 

0.15 

0.10 

0.09 

0.07 

30 

0.74 

1.49 

2.79 

4.89 

7.94 

0.015 

0.31 

0.24 

0.17 

0.15 

0.12 

40 

0.63 

1.27 

2.39 

4.19 

6.81 

0.020 

0.47 

0.37 

0.25 

0.22 

0.18 

50 

0.53 

1.06 

2.00 

3.49 

5.67 

0.025 

0.65 

0.50 

0.35 

0.30 

0.24 

60 

0.42 

0.85 

1.60 

2.79 

4.54 

0.030 

0.84 

0.66 

0.45 

0.39 

0.32 

70 

0.32 

0.64 

1.20 

2.09 

3.40 

0.035 

1.06 

0.83 

0.58 

0.49 

0.39 

80 

0.21 

0.43 

0.80 

1.40 

2.27 

0.040 

1.30 

1.02 

0.70 

0.61 

0.49 

90 

0.11 

0.21 

0.40 

0.70 

1.13 

0.045 

1.6 

1.22 

0.86 

0.72 

0.59 

100 

0 

0 

0 

0 

0 

0.050 

1.9 

1.4 

1.0 

0.84 

0.69 

0.06 

2.5 

1.9 

1.3 

1.12 

0.92 

0.07 

3.1 

2.4 

1.7 

1.4 

1.14 

0.08 

3.8 

3.0 

2.1 

1.8 

1.4 

0.09 

4.7 

3.7 

2.5 

2.2 

1.7 

0.10 

5.6 

4.3 

3.0 

2.5 

2.0 

0.11 

6.7 

5.1 

3.5 

3.0 

2.4 

0.12 

7.9 

6.1 

4.2 

3.5 

2.8 

0.13 

9.5 

7.2 

4.5 

4.2 

3.3 

0.14 

11 

8.4 

5.7 

4.8 

3.8 

0.15 

13 

10 

7.0 

5.6 

4.4 

0.16 

15 

11 

8.0 

6.6 

5.1 

0.17 

17 

13 

9.0 

7.5 

5.9 

0.18 

20 

15 

10 

8.4 

6.6 

0.19 

22 

17 

11 

9.7 

7.4 

0.20 

25 

19 

13 

11 

8.4 

0.22 

33 

25 

19 

15 

11 

0.24 

38 

29 

22 

17 

13 

0.26 

43 

33 

25 

20 

15 

0.28 

51 

39 

30 

23 

18 

0.30 

56 

43 

33 

25 

19 

0.32 

62 

47 

36 

28 

21 

0.34 

68 

52 

40 

31 

23 

0.36 

74 

57 

43 

34 

26 

0.38 

80 

61 

47 

36 

28 

0.40 

86 

66 

50 

39 

30 

\k 
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Step  1.  Reynolds  number  data  from  List  (1949)  along  with  F  factor 

data  from  Kinzer  and  Gunn  (1951)  were  inserted  in  the  expression 
1/2 

4irr(l  +  0.22  F  Re  )  to  extend  Table  C  2A  for  the  radius  interval  from 
0.21  cm  to  0.29  cm. 

Step  2.  From  radii  between  0.29  and  0.40  cm.  Re  was  calculated  by 
using  u,  the  temperature  dependent  kinematic  viscosity,  and  the  terminal 
velocity  values  given  in  Table  1.  The  F  factor  values  were  linearly 
extrapolated  from  the  data  given  in  Kinzer  and  Gunn  (1951) .  Entries 
in  Table  C.2A  were  then  computed  with  the  same  expression  as  given  in 
Step  1. 

b.  Discussion 

Kinzer  and  Gunn  (1951)  point  out  that  because  of  the  wide  range  of 
sizes  in  which  water  drops  are  found  in  nature,  the  physics  of  their 
evaporation  may  be  different  over  different  size  ranges.  They  suggest 
that  for  the  discussion  of  evaporation,  the  drops  should  be  treated  in 
three  following  categories: 

1.  Those  drops  that  lie  in  the  Stokes  Law  region  (r  <  0.004  cm); 

2.  Those  drops  in  the  radius  range  0.004  to  0.20  cm  that 

have  sufficient  speed  to  ventilate  adequately  the  transition 
layer  of  vapor  and  temperature. 

3.  Those  drops  with  a  radius  larger  than  0.20  cm  that  fall  so 
fast  as  to  be  deformed  and  whose  description  requires  special 
mathematical  analysis. 

Kinzer  and  Gunn  (1951)  worked  with  drops  in  the  second  category, 


and  this  research  extended  their  results  without  taking  drop  deforma- 
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tion  Into  account.  Watts  (1971)  calculated  mass  evaporation  rates  by 
using  heat  and  mass  transfer  analysis.  He  suggests  using  the  long 
radius  of  the  deformed  drop  rather  than  the  spherical  radius  in 
(C.  1).  This  suggestion  was  not  incorporated  in  this  study  because 
the  data  seemed  incomplete.  In  addition,  it  was  intended  to  compare 
the  results  obtained  in  this  study  with  those  of  other  workers  who 
used  the  original  Kinzer  and  Gunn  evaporation  data.  It  may  be  noted, 
however,  that  in  the  cases  treated  in  this  work,  collisional  breakup 
severely  limited  the  number  of  drops  with  radii  greater  than  0.2  cm. 


) 
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APPENDIX  D 

INVESTIGATION  OF  EVAPORATION  DATA 


a.  Background 

The  unexpected  Feature,  namely,  a  bimodality  in  the  drop  con¬ 
centration  at  the  lower  end  of  the  drop-size  range,  noted  in  Chapter 
5a,  is  examined  in  greater  detail.  Eq.  (53)  is  the  forward  differ¬ 
encing  formulation  used  to  calculate  drop  concentration  change  due 
to  evaporation.  For  all  initial  rainfall  rates  and  constant  downdraft 
velocities,  the  final  drop-size  distributions  were  similar  to  Fig. 

D.l.  The  bimodal  feature  in  these  distributions,  that  obviously 
is  due  to  evaporation,  is  not  documented  in  the  literature. 


E 

u 


r-> 


Radius  (cm) 


Fig.  D.l.  Bimodal  type  distribution  resulting  from  the 
evaporation  process. 

I 

i 
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In  order  to  Isolate  the  cause  of  this  trend.  Individual  tests  were 
performed  by  repeating  computations  In  which 

1.  The  number  of  drop-size  categories  was  decreased  from 
40  to  20; 

2.  Az  was  increased  from  25  to  50  m; 

3.  backward  as  well  as  forward  differencing  were  employed  in 
the  numerical  calculations ; 

4.  computation  was  performed  with  double  precision; 

5.  the  convergence  criteria  were  made  more  stringent  in  the 
iterative  scheme  for  obtaining  the  new  temperature  at 
each  level; 

6.  various  downdraft  velocities  (1,  5,  10,  and  15  m/sec) 
were  tried ;  and 

7.  various  initial  precipitation  rates  (25,  50,  75,  100 
mm/hr)  were  used. 

The  bimodal  charactertistic  appeared  in  the  final  drop  concentration 
for  each  of  these  computations.  The  final  drop  concentrations  exhibit¬ 
ed  a  decrease  in  the  radius  range  0.0053-0.0067  cm  and  increased  again 
to  a  maximum  in  the  0. 0187-0. 0191-cm  interval.  The  results  of  using 
wD  «  1  m/sec  showed  a  flattening  of  the  curve,  while  the  bimodality 
was  sharpened  by  stronger  downdrafts.  When  the  downdraft  was  held 
constant  and  the  initial  precipitation  rate  allowed  to  vary,  the 
bimodality  was  sharpened  by  larger  values  of  R.  When  the  number  of 
radius  categories  was  decreased  to  20,  the  decrease  in  final  drop 
concentration  occurred  in  the  radius  range  from  0.0088  to  0.0139  cm 
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with  a  peak  in  the  0.0173  to  0.0277- cm  radius  range.  No  difference 
was  noted  in  the  tests  in  2,  3,  4,  and  5  above. 

Following  the  above  tests f  the  sensitivity  of  the  Kinzer  and 
Gunn  (1951)  evaporation  data  was  investigated.  This  was  done  by 
smoothing  the  mass  evaporation  data  by  using  (C.  2)  as  reported  by 
Byers  (1965).  The  value  of  a’  and  b'  were  calculated  from  values 
of  the  various  parameters  at  the  mean  temperature  of  the  ralnshaft 
rather  than  the  temperature-dependent  values  used  initially.  This 
treatment  is  similar  to  that  used  by  Brazier-Smith  et  al.  (1973). 

Their  evaporation  results,  based  on  an  initial  Marshall-Palmer  distri¬ 
bution,  did  not  show  the  bimodal  characteristic. 

b.  Method  of  Computation 

The  following  values  were  treated  as  constant  in  the  subcloud 
region: 

T  -  285 

Ly*  2.47  x  106j/kg 

0.582  jm"1sec~1k"1 

D  ■  2.44  x  10  ^  m2/sec 

e  -  14.49  x  102  n/m2 
s 

p  *  1.445  x  10  5  m2/sec 

Table  117. B  in  the  Smithsonian  Meteorological  Tables  was  then  re¬ 
constructed  by  using  S-l(a'  +  b'_  with  a'  +  b'  *  4.003  x  102/sec-m, 
and  S  equal  to  the  saturation  values  in  Table  117. B. 

In  order  to  reconstruct  Table  117. A,  F  was  assumed  equal  to 


unity,  and  a  constant  value  of  u  was  used.  Terminal  velocity  was 
plotted  as  a  smooth  function  of  drop  radius  from  the  results  of 
Gunn  and  Kinzer  (1949). 

After  Tables  117. A  and  117. B  were  completed,  they  were  read 
into  the  evaporation  program.  The  results  did  not  exhibit  the  blmodal 
trend  depicted  in  Fig.  D.l,  but  rather  demonstrated  a  single  peak 
at  the  radius  of  0.0179  cm  for  w^  ■  10  m/sec,  and  at  0.0250  cm  for 
wD  ■  5  m/sec. 

c.  Discussion 

The  above  results  indicate  that  the  bimodality  in  Fig.  D.l 
is  inherently  caused  by  the  evaporation  data  of  Kinzer  and  Gunn 
(1951).  The  smooth  evaporation  results  are  consistent  with  the 
evaporation  results  of  Brazier-Smith  et  al.  (1973). 
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